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To Margarete and Ingrid 




Preface to the second edition 



Since the first edition of this book was published, already five years have passed, 
and during this period research on automotive control has flourished, and at 
the same time, the amount of industrial applications and products have also 
prospered. This means that there was a rich flora of possibilities regarding the 
selection of material when planning the second edition. In this process new topics 
have been added, important fields have been deepened, and in order to keep the 
number of pages down, sections of limited interest have been eliminated. 

In the selection of the new material we have concentrated on subjects that 
are both of current interest and importance, but at the same time are subjects 
that also contribute to a better understanding of basic processes and theories. 
The new material includes two completely new chapters: Diesel Engine Modeling 
and Diagnosis. 

In driveline control a new section on Anti-Jerk Control has been added. Large 
parts of Vehicle Dynamics and Control have been rewritten, which significantly 
improves the presentation of that material. Further, in this second edition, we 
have hopefully corrected most of the errors in the first edition, reviewed the 
nomenclature, and in order to facilitate to work with this book, added an index. 

The level of presentation has been thought trough to be suited for students at 
late undergraduate level or at early graduate level. The so called Bologna process 
is influencing universities in Europe at this moment, and in that perspective this 
book should be well suited for a course at the two year Masters level. 

We like to thank Dr. Dara Torkzadeh and Thomas Rainbow for their contri- 
butions in Diesel Engine Modeling, Dr. Mattias Nyberg for his work in Diagnosis, 
Julian Baumann for his participation in Anti-Jerking Control and Dr. Marcus 
Hiemer as well as Jorg Barrlro for their work in Vehicle Dynamics and reviewing 
this book. 

November 2004 



Uwe Kiencke 



Lars Nielsen 





Preface to the first edition 



Automotive control has become a driving factor in automotive innovation over 
the last twenty five years. In order to meet the enhanced requirements for lower 
fuel consumption, lower exhaust emissions, improved safety as well as comfort 
and convenience functions, automotive control had to be applied. 

In any area of technology, control design is an interplay between reality, 
physics, modeling, and design methods. This is also true in automotive control, 
and there has been extensive work done in research and development leading to 
a number of descriptions, models, and design methodologies suited for control. 



Goal of the book 

Our purpose of writing a book on Automotive Control is to present this interplay 
between thermodynamics, basics of engine operation, vehicle mechanics as well 
as parameter estimation and automotive control approaches. 

There are several good books available on the separate disciplines (some of 
the major references are in German). However, up until now there has not 
been a text available that explores more deeply the connections between reality, 
measurements, models and control design. 

It has been natural for us to treat all the major aspects of automotive control 
in the same book. This means that we cover engine, driveline, and complete 
vehicle. One reason is that there are similarities in methodology when analyzing 
and designing automotive control systems. This includes the point of view of 
finding models of suitable complexity and expressiveness. Another, perhaps more 
important, reason is that there is a strong trend that engine control, driveline 
control, and vehicle control rather than being separate will be more and more 
integrated, so that overall vehicle optimization is possible. 

It has also been important to us to show real measurements. This gives 
a reader the possibility to see how models are approximations of reality, and 
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to judge the modeling assumptions. A consequence of this approach is that 
we have selected to treat systems that are close to some of those utilized in 
actual vehicles, rather than discussing speculative systems or presenting purely 
theoretical results. 

Intended readers 

This book should enable control engineers to understand engine and vehicle mod- 
els necessary for controller design and should introduce mechanical engineers into 
vehicle-specific signal processing and automatic control. 

In fact, our inspiration to write the book came from this. We are both mem- 
bers of the IFAC technical committee on Automotive Control (with the first 
author being the chairman). We met there and also at SAE meetings, and we 
saw the potential value of bridging a gap that was obvious to us. However, even 
more important to us is to share some of the fun and excitement that goes into 
the area of Automotive Control Systems and thus give it the attention it deserves. 

Organization of the book 

The outline of the book starts with engines, continues with drivelines, and finally 
deals with the vehicle. 

Chapters 2 to 4 treat engines with regard to basics, thermodynamics, models, 
control, and advanced concepts. All the major control systems and their design 
are treated. The thermodynamic models in Chapter 2 deal with parameters that 
vary under one cycle and the resolution of interest is typically one crank angle 
degree, whereas the time scales of mean value models are in the order of 1 to 
several engine cycles, and the variation in variables that are considered are also 
averaged over one or several cycles. These models form the basis for understand- 
ing the complex phenomena that influence the engine operation, efficiency and 
emissions. They also serve the purpose of describing the properties influencing 
control design and performance in Chapters 3 and 5. 

The driveline (engine, clutch, transmission, shafts, and wheels) which is a 
fundamental part of a vehicle is the topic in Chapter 7. Since the parts are elastic, 
mechanical resonances may occur. The handling of such resonances is basic 
for functionality and driveability, but is also important for reducing mechanical 
stress and noise. Two important modes of driveline control that are treated are 
driveline speed control and driveline torque control, having their applications in 
cruise control and automatic gear shifting control. 

Vehicle dynamics control systems help the driver to perform the task of keep- 
ing the vehicle on the road in a safe manner. These systems are thus often 
safety-oriented, which means that they only interact in situations where they 
can reduce the possibility of an accident, but then they affect the immediate 
behavior of the vehicle within fractions of a second. Some systems are also used 
for improving the comfort of the driver. The performance of a vehicle, regarding 
the motions coming from accelerating, braking, cornering, or ride, is mainly a 
response to the forces imposed on the vehicle from the tire-road contact. Much of 
study of vehicle dynamics is a study on why and how these forces are produced 
and how they can be effectively understood and treated in simplified models. 
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The basics of these models and some associated control systems are presented in 
Chapters 8 to 10. 

Chapter 11 is the exception from that all the systems and principles in this 
book is close to some of those utilized in actual vehicles. The reason is that road 
and driver modeling is part of simulation design rather than part of a vehicle. 
Nevertheless, it is important to realize that road and driver models are important 
parts in the design cycle of automotive systems design due to the importance of 
advanced simulation. 

Background and use of the book 

The material in this book has been used in courses at the universities of Karls- 
ruhe, Germany and Linkoping, Sweden. It is well suited for the later stages 
(third or fourth year) of the engineering programs at our technical institutes 
( “Diploma-engineer” , “Master of Science” ) . 

The book, to a large extent, covers the basic material needed, but of course 
it is advantageous to have a background from basic undergraduate courses in 
automatic control, signals and systems, mechanics, and physics. 

The course lay-out includes problem-solving sessions and laboratory exper- 
iments. The laboratory assignments typically include measurements, building 
models of the type treated in the book, and finally designing controllers and 
simulating them. Here students with more background, for example in modern 
control, can do more elaborate designs. This is also the case when the book is 
used in an introductory graduate course. 

The authors 

Dr. Kiencke’s experience in this field started in the early nineteen seventies 
when developing adaptive lambda control and knock control at Robert Bosch 
Corporation. In the following years more complex approaches for engine mod- 
elling [2], [22] and controller design [63] were published. At that time he headed 
a team that developed the vehicle communication network ’’Controller Area Net- 
work (CAN)” [67]. Networking allowed to combine formerly stand-alone control 
schemes into an integrated vehicle control system. In the early nineteen nineties 
he joined the University of Karlsruhe in Germany where he could intensify engine 
and vehicle dynmaic control research. 

Dr. Nielsen has more than fifteen years background in academic mechatronics 
research (obtaining a good start at the Department of Automatic Control in 
Lund, Sweden) . He has during that time continuously collaborated with industry, 
and has lead joint research projects with Scania AB, Mecel AB, Saab Automobile 
AB, Volvo AB, and Daimler Chrysler. He is since 1992 holder of the chair Sten 
Gustafsson professor of vehicular systems at Linkoping University in Sweden. 

Acknowledgments 

The control systems presented were mostly developed within a team. Therefore 
the first author would like to thank especially the following cooperation partners: 
Dr. Martin Zechall in lambda (air-fuel ratio) control, Dr. Boning in knock control 
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and engine map optimization, Alfred Schutz in engine idle speed control, Heinz 
Leiber in ABS braking control, Dr. Michael Henn in misfire detection, Dr. Achirn 
Daiss in vehicle modelling and identification and Dr. Rajjid Majjad in road and 
driver modelling. It was a great pleasure to cooperate with these people and it 
created many friendships. The second author is especially indebted to Magnus 
Petterssson for joint work in driveline control, and to Lars Eriksson for joint work 
in engine modelling and control. Also Lars-Gunnar Hedstrom, Jan Nytornt, and 
Jan Dellrud deserves special mentioning as research dedicated industrial partners. 

Furthermore we both thank Christopher Riegel, Jochen Schontaler, Dara 
Torkzadeh, and Dr. Tracy Dalton for their tremendous effort to translate and 
revise parts of the book, as well as Dr. Dietrich Merkle as a publisher. 

Last but not least we to thank our families and especially our wives Margarete 
and Ingrid for tolerating that so much weekend and vacation time was dedicated 
to this book. 



Being in November 1999 looking forward to the next millennium, we hope that 
readers will share some of the excitement that comes along with Automotive 
Control Systems. 



Uwe Kiencke 



Lars Nielsen 
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1 Introduction 



Vehicles are now computerized machines. This fact has had an enormous effect 
on the possibilities for functionality of vehicles, which together with needs and 
requirements from customers and from society have created vigorous activities in 
development. 



1.1 Overall demands 

The overall demands on a vehicle are that it should provide safe and comfort- 
able transportation together with good environmental protection and good fuel 
economy. This means that there are three main objectives for automotive control 
systems: 

• Efficiency, which leads to lower fuel consumption. 

• Emissions should be low to protect the environment. 

• Safety is of course a key issue. 

There are a number of additional objectives like comfort, driveability, low wear, 
availability, and long term functionality. 



1.2 Historic remark 

Many of the technologies that today are considered advanced, sometimes even 
new, have been around for a long time. It is therefore interesting to ask our- 
selves why these technologies are surfacing now as commercial products. Direct 
injection of gasoline engines is one example. These concepts are not new, even 
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if sometimes presented so, but the novelty is instead that they now with proper 
control can achieve competitive functionality and performance. 

It is thus the breakthrough of computer control that is a driving factor. A 
good example is ABS (Anti-lock Braking Systems) which is an old idea, but it 
was not functional enough using mechanical solutions or analog electronics. Now 
these systems are readily available and widely spread. 

1.3 Perspectives 

Looking at the future, the three overall objectives above will be in focus. The 
demands of reduced emissions and advanced diagnosis functionality are steadily 
increased by legislators and customers. The key areas that help meeting the 
increased demands are the development of control and diagnosis functions in 
the control units. Further, this functionality will have to be obtained not only 
when the car is new, but over a sustained period of time. Other examples are 
improved stability due to handling control, and improved driveability due to 
driveline torque control, which also can be used to e.g. reduce clutch wear. 

Regarding methodology development, mathematical models will play an im- 
portant role. They will be used for model based control and diagnosis. They will 
also be the basis for e.g. sensor fusion, adaptive control, and supervision. 

Co-design 

Automotive control will not only improve existing vehicle designs. It will also to 
a large extent change the view of vehicular systems design leading to: 

• New mechanical designs. These new designs are made possible by, and rely 
on, the existence of a control system. 

Design of vehicles is thus evolving into co-design of mechanics and control. The 
goals for this development can be set high, and the perspectives on automotive 
control systems are therefore concluded with an inspiring mind teaser: 

A mind teaser 

It can not be ruled out that a car can function as an air cleaner for usual town air. 
In a typical town in the industrial world, the air is typically somewhat polluted 
from many sources, including cars, but also due to house heating and industries. 
Existing and upcoming technology lowering exhaust emissions are such that the 
concentrations in the exhaust after the catalyst can be lower than in town air. 
This means that the originally available pollutants have been combusted or have 
been collected in the catalyst. 

Since such a perfect combustion produces only water and carbon dioxide, the 
problem of carbon dioxide is then a matter of less fuel consumption. 

If at all possible, such a development will rely heavily on automotive con- 
trol systems since the car has to function as an air cleaner under all possible 
circumstances, load variations, and driving styles. This is a mind teaser, but 
understanding the material in this book is a first step. 




2 Thermodynamic Engine Cycles 



In this chapter, the thermodynamic characteristics of basic engine cycles are ex- 
plained. For each concept, the thermal efficiency is derived from thermodynamic 
equations. 



2.1 Introduction to Thermodynamics 

A short introduction to thermodynamics is given in this section. It is essential to 
understand the thermodynamic processes for the different types of engine models 
explained in this book. 

2.1.1 First Thermodynamic Law 

Ideal gases are always in a gaseous state and they behave according to the ideal 
gas equation : 

pF = ruRd (2.1) 

where: 

p is the pressure in N/m 2 
V is the volume in m 3 
m is the mass in kg 

R is the gas constant R = 287.4 m 2 /(s 2 K) 
d is the absolute temperature in K 

The expression p V has the unit Joule and it depends on the mass and temper- 
ature of the gas. 

The thermal energy is defined as: 



q = cm d 



(2.2) 
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where: 

q is the thermal energy in J 
c is the specific heat constant in J/{kg K) 



This equation is valid for all kind of materials in solid, liquid and gaseous state. 
The specific heat constant is material depending. In gases, different state changes 
have to be distinguished which will be discussed in the next section. 



First Thermodynamic Law 



dq = du + dw (2.3) 

dq is the differential change of thermal energy. Energy which is 
brought into the gas is positive. 
du is the differential change of internal energy. Internal energy 
which is brought into the gas is positive. 
dw is the differential change of mechanical work. Mechanical work 
which is brought into the gas is negative. 

Please note, that this may be defined differently in other books. 

The thermal energy dq brought into a gas of constant mass results in a rise in 
it’s internal energy du and/or is transformed into mechanical work outside the 
gas. Thus, energy is neither created nor destroyed. This equation is solved by 
integration to compute the energies in state changes. 



Volume Change 

There are two basic principles how mechanical work can be delivered to or from 
the system: By a change of volume or a change of pressure. 

In engines, a limited amount of compressed gas expands in a cylinder. The 
differential mechanical work dw v depends on the force F upon the piston and 
the differential piston stroke ds: 

dw v = F ds = p A ds 



where A is the cross-section of the cylinder and p the pressure. More general, 
this yields to: 

dw v = p dV (2-4) 

The work dw v is equivalent to the kinetic energy brought into the mechanical 
system due to the expansion of the gas in the cylinder. Mechanical energy brought 
into the gas is negative: dV < 0 (compression) and therefore dw v < 0. Output 
of mechanical work from the gas is positive, dV > 0 (expansion), dw v > 0. This 
can also be seen in Figure 2.1. 
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Figure 2.1 Mechanical work due to volume change. 




Figure 2.2 Mechanical work due to pressure change 



Pressure Change 

In turbines, kinetic energy is transmitted by a continuous gas flow between two 
locations of different pressure: 



dw p = —V dp (2-5) 

The work dw p is equivalent to the loss of potential energy of the gas. Figure 2.2 
graphically explains the state change in the pV-diagram. Output of mechanical 
work from the gas is positive because the pressure drop is negative: dp < 0 
(expansion), dw p > 0. 

The two different types of work are graphically explained in Figure 2.3 and 
they are linked by the following equation: 
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Figure 2.3 Relationship between the two mechanical works. 



2 2 

W v = J p dV = - J V dp - PiVi + p 2 V 2 (2.6) 

l l 

Assuming a constant mass and a constant temperature, the two energies are 
equivalent as piV\ = P 2 V 2 according to the ideal gas equation. 

Enthalpy 

The state of the gas can be described with several state variables: 
u Internal energy of the gas. It is equivalent to the amount 

of thermodynamic energy due to the movement of the atoms 
at a given temperature. 

p V State of the gas. The state depends on the volume and pressure 
of the gas. It is related to the internal energy of the gas. 

A state change can be caused by a change of pressure, a change of volume or 
both: 

d{p V) = p dV + V dp (2.7) 

p dV Kinetic energy caused by a change of volume. 

V dp Potential energy caused by a change of pressure. 

The enthalpy is another state variable of the gas: 

h = u+p V (2.8) 

The differential is: 

dh = du + p dV + V dp (2.9) 

Compared to the first thermodynamic law in Equation 2.3 



dq = du + dw 
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the thermal energy can be written as: 

dq = dh — V dp = du + p dV (2.10) 

For isobaric expansion, we have dw v = p dV and V dp = 0 and therefore dq = dh. 
This yields to: 

dq = du + pdV (2-11) 

For isochoric pressure drop, we have dw p = —V dp and p dV = 0 and therefore 
dq = du. This results in: 

dq = dh~Vdp (2-12) 

The enthalpy h. can be very useful as a state variable when dealing with gas 
turbine engines. The internal energy u is useful when dealing with piston engines. 



2.1.2 Specific Heat Constant 

Input of Thermal Energy at Constant Volume 

A temperature rise can be observed when increasing the thermal energy q of a 
constant mass of gas at a constant volume. The specific heat constant c v is 
defined as: 



1 



Cv — 

m 



The differential of the internal energy is 



du — 



&u\ 

w) t 



dV 



(du 

\dd 



dd 



Inserting this into the first thermodynamic law 2.11yields: 



f du 



\dV 

(du 



(du 



d( l = ( 7iT7 ) dv +\ vrj dd+pdV 



\dd 



\dd 

(s' * 



= Us *■&+ hwr +P UV 



du 



dV 



At a constant volume we have dV = 0. After division by dd: 

dq\ / du 



ddj 



v 



dd 



Equation 2.13 can be written as: 



1 / du 
m \dd 



(2-13) 



(2.14) 

(2.15) 
(2-16) 



(2.17) 



We have shown that a change in the thermal energy dq at constant volume is 
equivalent to a change in the internal energy du of the gas. 
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Input of Thermal Energy at Constant Pressure 



A temperature rise can also be observed when increasing the thermal energy of 
a constant mass of gas at a constant pressure. Please note the difference to 
the previous section, as the specific heat constant c p is now defined as: 



1 



dq 

dd 



(2.18) 



The differential of the enthalpy is: 

MS 



d0+(^ 

dp 



dp 



This relationship can be inserted into 2.12 and then be simplified because dp = 0: 



, ,'dh 

d d = M 



dd 

fdh 

\dd 



dd + 



dd + 



dp-V dp 

d PJn 
fdh, 

\dp 



— V ) dp 



- 'S' « 



The last equation is divided by dd: 

dq 
dd 



dh 

dd 



Similar to the above section, Equation 2.18 can also be written as: 



1 fdh 



Cr, — 



m \ dd 



(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 



( 2 . 22 ) 



The enthalpy of the system is changed. The input of thermal energy dq at 
constant pressure not only increases the internal energy du, but also produces an 
output of work p dV. 



Gas Constant and Adiabatic Coefficient 

The relationship between c v and c v can be seen by using the enthalpy (Equa- 
tion 2.8) and the ideal gas equation (Equation 2.1): 



h 


= u + pV 


(2.23) 


dh 


= du + d(pV) 


(2.24) 


in c p dd 


= m c v dd + mR dd 


(2.25) 


Cp 


= Cy + R 


(2.26) 



The ideal gas constant R is equal to the difference between c p and c v : 

R = c p — c v 



(2.27) 
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Figure 2.4 Isothermic state change 



The ratio of c v to c v is called the adiabatic exponent: 




Cy 

Under normal conditions its value for air is n a ir = yf = 1-4. 



(2.28) 



2.1.3 State Changes of Ideal Gases 

Isothermal Change: 0 = const 

Assuming a constant temperature, the ideal gas equation pV = m Rd yields: 

rn R d const 



where p V = const. The state change is equivalent to a hyperbolic curve in the 
pU-diagram which can be seen in Figure 2.4. The isothermal state change is 
characterized by using the ideal gas equation: 

PiV 1 =p 2 V 2 (2.29) 



The slope of the curve in the pU-diagram is: 

dp p 

dV = ~V 

The internal energy of the gas remains constant: 



(2.30) 



du = mc v dd = 0 



Hence, the first thermodynamic law can be written as: 



dq = du + dw = dw 
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The enthalpy also remains constant because of p V = const and d{p V) = 0: 

dh = du + d(p V) = 0 

The kinetic energy caused by the state change can be derived by using: 

d{p V) = p dV + V dp = 0 



Hence, it is simplified to: 



and integrated: 



dw = p dV = —V dp 
2 2 

tt'i 2 = [ p dV = — [ V dp 



Finally, by using the ideal gas equation: 



! 1 V 2 

u>i,2 = Qi, 2 = m R d / — dV = m R d In — 



V 



Vi 



Work w 12 is generated from the system by expanding the gas. The same amount 
of thermal energy must be brought into the system in order to compensate 
for the work and to maintain the constant temperature. 



Isobaric Change: p = const 

A state change at constant pressure is described by the ideal gas equation: 



(2.31) 



mRd 

p = — — — = const 

and therefore: 

0i = Vi 

d 2 V 2 

The state change is equivalent to a horizontal line in the pH-diagram which can 
be seen in Figure 2.5. The isobaric state change 1 — > 2 results in an output of 
mechanical work Wi t2 due to the expansion of the gas, and at the same time it 
results in an increased internal energy u due to the temperature rise. The input 
of thermal energy q is equivalent to the increase of the enthalpy h. 

In this case, the first thermodynamic law can be written as: 



dq = dh — V dp = dh 



since pressure change vanishes: dp = 0. The differential change of thermal energy 
is given by: 

dq = to c p dd 



Integration leads to: 



A 

<7i,2 = to c v J dd = to c p (d 2 — di) 



(2.32) 
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Figure 2.5 Isobaric state change 



The mechanical work can be derived by using the fact that dq = dh and dh = 
du + p dV : 

dq = du + p dV 

Because of dw v = p dV , the change of kinetic energy is: 



p dV = dh — du 


(2.33) 


TO Cp did — TO c v did 


(2.34) 


mR did 


(2.35) 



Integration leads to: 



Z 

wi t 2 = ruRD J did = mR(id 2 — $ 1 ) 



Please note that only a small portion of thermal energy q is converted into kinetic 
energy w. The ratio of the input of thermal energy to the output of mechanical 
work can be expressed: 



dw m(c p — c v ) did c p — c v k — 1 

dq mc p did c p k 



which is 0.29 for air. 



Isochoric Change: V = const 
In this case, the ideal gas equation is: 



V = 



mRid 

P 



= const 
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Figure 2.6 Isochoric state change 



and therefore: 



$1 _ Pi 
d 2 P2 



(2.36) 



The state change is equivalent to a vertical line in the pV-diagram illustrated in 
Figure 2.6. The mechanical work vanishes for an isochoric state change 1 — > 2. 
The input of thermal energy dq is equivalent to the increased internal energy du. 



As dq = du + p dV = du, the thermal energy is given by: 



dq = mc v dd 
2 


(2.37) 


51,2 = rn c v f dd = mc v {d 2 — d\) 


(2.38) 



l 



The enthalpy is increased more than the internal energy: 

dh = du + p dV + V dp 

where p dV = 0 and V dp = m R dd. The change of enthalpy is therefore: 
dh = m c v dd + mR dd — m c p dd 

The internal energy is du = dh — V dp. It has the same value as the thermal 
energy: 

du = m(c p — R) dd = m c v dd = dq 



Isentropic or Adiabatic Change: q = const 

Assuming an insulated gas volume, the thermal energy of the gas remains con- 
stant. In an isentropic state change 1 — > 2, the gas expands and mechanical work 
is delivered from the gas while its internal energy is reduced. The gradient of the 
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Figure 2.7 Isentropic state change 



adiabatic change in the pH-diagram is steeper than the hyperbolic curve of the 
isothermic change. This is due to the fact that no thermal energy is brought into 
the gas. 

dq = du + p dV = 0 

Hence, the first thermodynamic law can be written as: 



— du = p dV 



This means that the mechanical work is equivalent to the loss of internal energy. 
The change of enthalpy is 

dh = du + p dV + V dp = V dp 

because of dq = du + p dV = 0. The isentropic change can be considered as the 
sequence of two state changes: 

• An isobaric change where dq p = m c p dd p = dh. 

• An isochoric change where dq v = mc v dd v = du. 

The summation of exchanged thermal energies is supposed to be zero: 

dq = dq p + dq v = 0 (2.39) 

c p dd p + c v dfl v = 0 (2.40) 

From the isobaric and the isochoric state changes we get: 

isobaric : dd p = 



p dV 
mR 
V dp 



(2.41) 

(2.42) 



isochoric : did. 



mR 
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This can be inserted into Equation 2.40 



p dV 
Cp mR 



+ c v 



V dp 
to R 



= 0 



which can be rearranged to obtain the gradient in the pV-diagram: 



dp _ _CpP^ _ _ p 

dV~ c v V~ V 



(2.43) 



The differential pressure change for isentropic state change is increased by the 
factor of k compared to the isothermic state change (see Equation 2.30). 
Equation 2.43 can be integrated: 



hip = —k In V + const 



This yields: 

pV K = const 

An isentropic change 1 — > 2 is characterized by: 

PiVi = p 2 V 2 

The kinetic energy 'fCi j2 is obtained by integration: 



(2.44) 



Z 

wi,2 = j P dV 



Exchanging p with p = const V K leads to: 

2 



W! 2 = const [ V~ K dV = - V^~ K ) 

J 1- K 



Replacing const = p{V{". 



wi,2=PiVf-^—(V 2 1 - K - V 1 1 ~ K )=p 1 V 1 -^— (7 

1 — K 1 — K V \ Vi ) 



- -1 



Further simplifications: 



(Vi\ K 1 = P 1 VIP 2 V 2 _ P 2 V 2 _ ■& 2 
V^2/ P2V 2 k piVi p{V 1 ill 

The work is given by: 



1 

1 — K 




w 1,2 = Pi El 




2.1. INTRODUCTION TO THERMODYNAMICS 



15 



By inserting the ideal gas equation p\V\ = mRDi = m(c p — c„)t?i and the fact 
that k = Cp/c v , which yields: 

w\,2 = m(cp - c v )D i — ^ ^ ^1 - (2.45) 

= m(c p - Cy) c„ \ (#i ~ # 2 ) (2.46) 

C v 

and finally, the mechanical work for an isentropic state change is given by: 

101,2 = mc„(i?i — $ 2 ) (2.47) 



Polytropic Change 

Polytropic state changes are necessary to model a realistic isentropic process with 
insufficient insulation. Therefore, the polytropic exponent n has values smaller 
than the isentropic exponent k: 



p V n = const 



with n < k. 

2.1.4 Thermodynamic Cycles 

Entropy 

Two different types of state changes have to be distinguished: reversible and 
irreversible state changes. They can be described by the following characteristics. 
Reversible state changes: 

• enable the system to return to the initial state 

• do not need energy from outside the system when restored to their initial 
state. 

• do not leave a permanent state change in a closed system after restoration. 
Irreversible state changes: 

• go into one direction. For example, expansion of molecules when a larger 
volume is available to them 

• need energy from outside which is transformed into thermal energy. For 
example, a falling stone. 

• do not return to the initial state. 

The entropy is defined as: 

S = l (2-48) 

An example is given to explain the entropy: 
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Figure 2.8 Ideal cycle process 



There are two bodies at two different temperatures i9. Thermal energy q is 
exchanged by conduction from the hot body at temperature to the cold body 
at temperature 1 I 2 • The change of entropy within the system after the transfer 
of thermal energy can be described by: 

AS =f~f (2 - 49) 

Even when qi=q 2 there is a positive change of entropy AS due to the irreversible 
process of energy flow from the hot to the cold body. Thermal energy can only 
be transferred from hot to cold media and not in the opposite direction. The 
entropy of the cold body 2 is larger than the entropy of the hot body 1 as the 
entropy is inversely proportional to the temperature. 

A reversible state change is characterized by: 

AS = 0 (2.50) 

and an irreversible state change is characterized by: 

AS > 0 (2.51) 

which is also called the second thermodynamic law. It means that the entropy 
within the system always increases for real state changes. 



Ideal Cycle Process 

Assuming reversible state changes, the following ideal cycle process can be de- 
scribed: The gas absorbs thermal energy qi in an isothermal and reversible state 
change at high temperature $ 1 . The gas then loses the thermal energy (72 in an 
isothermal and reversible state change at low temperature '&2 ■ Kinetic energy 
w 1,2 is delivered by such an ideal engine which is equivalent to the difference of 
the thermal energies: W \2 = q\ — < 72 - The entropy change shall be zero: 



92 _ 9i 
$2 1? 



and therefore: 



9i 92 
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Figure 2.9 pE-diagram and dS-diagram of Carnot Cycle. 



The thermal efficiency of this ideal cycle process is equivalent to the ratio of 
mechanical work w 1,2 to the absorbed thermal energy q-\ : 



wi, 2 _ gi - ff2 _ $1 - # 2 _ 1 _ $2 
q 1 q 1 $1 ■di 



(2.53) 



It can be seen that the thermal efficiency 77 only depends on the ratio of the 
absolute temperatures di to $2- ?? is always smaller than 1. Please note that in 
non-ideal cycle processes the entropy change is always non-zero, i. e. AS > 0. 



Carnot Cycle 

The Carnot cycle is characterized by four state changes. At the end of the process, 
the initial state is reached again which is the case for all periodic processes in 
engines. Two different cycles have to be distinguished: 

Reaction within a closed combustion chamber: Engines with periodic 
combustion cycles and emission of mechanical work. 



Gas flow through an open combustion chamber: Turbines with continu- 
ous combustion and emission of mechanical work. 



The Carnot cycle is used as reference model for thermodynamic processes. It is 
illustrated in the pE-diagram of Figure 2.9: 

1 — > 2 Isothermal expansion of the gas from El to V2 while the temperature 
remains constant di = $2- 

V 2 

wi j2 = Vi,2 =mR $ 1 l n 7 j~ 

VI 

The mechanical work u>i ; 2 is delivered since the gas expands. Thermal 
energy gi^ of the same amount has to be brought into the system in order 
to keep the temperature constant. 



2 — > 3 Isentropic expansion of the gas from V 2 to V3. No thermal energy is 
exchanged g2,3 = 0. Output of kinetic energy: 

W2,3 = m c v (i?2 - d 3 ) 
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3 — y 4 Isothermal compression of the gas from V 3 to V 4 while the temperature 

remains constant $3 = ^ 4 . The mechanical work 

V 4 

W3,4 = 93,4 = mRv 3 In — 

*3 

is needed which is negative in this case because it is delivered into the gas 
(V4 <V 3 ). 

4 — > 1 Isentropic compression of the gas from V 4 to V\. No thermal energy is 

exchanged 92,3 = 0. The mechanical work is also negative in this case 

($1 > tf 4 ): 

u>4,i = mc v {D 4 — $ 1 ) 

Now, all kinetic energies of the Carnot cycle have to be added to derive the 
thermodynamic efficiency: 

W = Wl,2 + w 2) 3 + ^3,4 + W4,l (2.54) 

= mRD 1 In — ^ + mc v (d 2 — $ 3 ) + mRi) 3 In + mc,(i ) 4 — i?i)(2.55) 
Vi V 3 

The work u> 2,3 and w 4) i of the isentropic state changes compensate each other 
because of di = $ 2 and $3 = i? 4 : 

Vo Va 

w = m R(il 1 In — + $ 3 In — ) 

Vi I /3 

To simplify this equation, the characteristic equations of the isentropic process 
1 — > 2 and 3 — > 4 are used: 



P V K = P V-V K ~ 1 


= mRD-V K 1 = const 


VzV*- 1 = 1?3^3 K_1 


, ^y*- 1 = ViVf- 1 


my - 1 


V 1 


w 


t?4 



(2.56) 

(2.57) 

(2.58) 



Because of $ 4 = $2 and $3 



and therefore: 



= D 4 we get 



55 

?v 

1 

II 


my 1 

m; 


(2.59) 


V 2 

V) 


V 3 

v 4 


(2.60) 


ln ^4 = 


" ln ^ 


(2.61) 


w = m 7?(-$i 


K, 

-^ 3 )iny 





This result is now used for the thermodynamic efficiency: 



w _mR{tf i-tf 3 )ln^| 
91,2 mRH iln|| 

1 - ^3 _ , $3 



(2.62) 



(2.63) 
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The efficiency depends only on the temperature ratio of $3 to $ 1 . 
compression ratio £ : 



£ = 



v± 

b 



Using the 



(2.64) 



the efficiency is ($3 = $ 4 ): 



1 



1 

1 



04 

01 

1 



1 



£ k ~ 



1 




(2.65) 

( 2 . 66 ) 



The thermal efficiency can be explained graphically in the pU-diagram in Fig- 
ure 2.9: It is equivalent to the ratio of the integral within the cycle to the area 
which is produced by integrating qi^- The efficiency is improved by increasing 
the area within the cycle. This can only be achieved by a rise of the compression 
ratio e or the temperature ratio. 

To give an example: Assuming absolute temperatures: 



(2.67) 

( 2 . 68 ) 

the thermal efficiency would be: 



= 2800 ° K 
i9 3 = 300 °K 



77 = 0.89 



and the compression ratio: 

£ = 266 

supposing a k = 1.4. Such a compression ratio can hardly be produced by real 
engines. Typical compression ratios are at least ten times smaller. 



2.2 Ideal Combustion Engines 

Commonly used combustion engines in cars are four-stroke engines. They have 
two intermittent cycles: the gas is compressed, combusted and expanded in the 
first cycle, and the gas is exchanged in the second cycle. In this section the second 
(or passive) cycle will not be considered to simplify the mathematical derivations. 
The processes related to the second cycle will be discussed in Chapter 3. 

Two different types of combustion engines have to be distinguished: 

1. Spark- ignited (SI) Engine: Combustion caused by an electric spark- ignition. 

2. Diesel Engine: Combustion caused by self inflammation due to compres- 
sional heat. 

In most sections, p represents the in-cylinder pressure, V the cylinder volume, il 
the in-cylinder temperature, S the entropy, q the thermal energy of the gas, u its 
internal energy and h its enthalpy. 
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Figure 2.10 pF-diagram (left) and d.S'-diagrarn (right) of the SI engine process 



2.2.1 Spark-ignited (SI) Engine 

The first SI engine was presented by Nikolaus Otto in 1862. The combustion 
process can be modeled as an isochoric process where the gas volume is con- 
sidered to be constant. The pF-diagram in Figure 2.10 illustrates that the gas 
volume does not change between step 2 and step 3. The ratio of maximum to 
minimum volume is given by: 



e = 



v 2 



(2.69) 



This ratio e is called the compression ratio of the engine. The different steps 
for a complete cycle in the pV - diagram and in the dS-diagram can be seen in 
Figure 2.10. Mathematically they can be described as followed: 



1 — > 2 : Isentropic compression, dq = 0: 



dq 

91,2 

dw 

U>1,2 



du + dw = 0 

0 



— du = — m c v dd 
2 

— J mc v dd = —to c v (d 2 






The work Wi t 2 is used to compress the gas and therefore, it is negative. 
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2 — > 3 : Isochoric input of thermal energy, dV = 0: 

dw = p dV = 0 

3 

u > 2 ,3 = J pdV = 0 
2 

dq = dii = m c v dd 
3 



92,3 = mc v J dd = mc v (d 3 - d 2 ) 

2 

The increased thermal energy g 2 , 3 is caused by combustion of the gas. 



3 — > 4 : Isentropic expansion, dq = 0: 



93,4 

dw 



= 0 



— du = —mc v dd 

4 



W 3 ,4 = 



— J m c v did = —m c v ($4 — $3) 



This state change describes the power stroke of the engine where is 
the output of kinetic energy from the gas, which is positive ($4 < i? 3 ). 



4 — > 1 : Isochoric heat loss, dV = 0: 

dw = p dV = 0 



«> 4,1 

dq 

94,1 



pdV = 0 



du + dw = mCyd'd 
1 

mc v J dd = mc v ( i?i 



0 4 ) 



4 

The loss of thermal energy q^i is due to the gas exchange: The burnt hot 
gas is pumped into the exhaust and the combustion chamber is filled with 
a cold mixture of unburnt fuel vapor and air (q^i is negative because of 
$1 < $4). 



The thermal efficiency of the engine is equivalent to the ratio of all the kinetic 
energies to the input of thermal energy g 2 3 at the combustion of a complete 




22 



2. THERMODYNAMIC ENGINE CYCLES 



cycle: 



Vth 



w 1,2 + w 2 , 3 + UJ 3A + W 4; 1 

92,3 

mc v (—'i 9 2 + $i - $4 + $3) 
mc v (fl 3 - $ 2 ) 

$4 - i?i 
^ 3-^2 
^1 ^4/^1 - 1 
$2 ^3/^2 - 1 



The relationship for isentropic changes 1 
the equation: 

t ± = (VsY ~ 1 

$ 3 v^J 



2 and 3 — > 4 can be used to simplify 



J_ = l 1 

e K_1 ^2 



(2.70) 



This yields: 

Vth = l-^ (2.71) 

Please note that the thermal efficiency rjth does not depend on the absolute 
temperature values. It mainly depends on the compression ratio e. Example: 
For a compression ratio of £ = 11 and an adiabatic coefficient of k = 1.4 the 
theoretical thermal efficiency pth is: 



Vth = 0.617 



2.2.2 Diesel Engine 

Rudolf Diesel developed this engine from 1893 to 1897. In a diesel engine, the 
combustion takes place in an isobaric state change during the downward move- 
ment of the piston. At the beginning of this process the combustion is controlled 
by the injection of fuel to maintain a constant pressure at the expansion from 2 to 
3. The isobaric state change is indicated between steps 2 and 3 in the pE-diagram 
in Figure 2.11. The more fuel is injected, the longer the distance between steps 2 
and 3 and the larger the volume ratio: 



V 3 = #3 
V 2 ‘d 2 



(2.72) 



This ratio is called injection ratio or load. The injection ratio p has an impact 
on the thermodynamic efficiency which is derived after explaining the different 
parts of the cycle: 



1 — ■> 2 : Isentropic compression, dq = 0: 



dq = du + dw = 0 
91,2 = 0 

dw = —du=~mc v d r d 
wi ,2 = -mc v {’d 2 -'di) 
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q 2.3 




Figure 2.11 pV-diagram for Diesel Engine 



The mechanical work w\^ is used to compress the gas (equivalent to the SI 
engine). It is negative. 

2 — > 3 : Isobaric gain of thermal energy, dp = 0: 

dq = dh — V dp = to c p did 
< 72,3 = mc p (’d 3 -i? 2 ) 

dw = p dV = mRdd 
W2,3 = m R(d3 - d 2 ) 

In this process, the combustion generates the thermal energy ( 72,3 and pro- 
duces the kinetic energy u> 2 , 3 - 

3 — > 4 : Isentropic expansion, dq = 0: 

<73,4 = 0 

dw = —du=—mc v d , d 

w 3t 4 = —mc v (il4 — $ 3 ) 



Note that 103,4 is positive since $4 < $ 3 . 
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4 — » 1 : Isochoric heat loss, dV = 0: 



dw = p dV = 0 

l 

u>4.i = J pdV = 0 

4 

dq = du + dw = mc v dd 

l 

<Z4.i = mc v J dd — mc v (di — $4) 



Note that 94,1 is negative since d\ < $4. 

With k = — and R = (c p — c„) the thermodynamic efficiency of the diesel engine 
can now be calculated: 



Pth = 



101,2 + U>2,3 + U>3A + W4,l 
92,3 

-mc v (d 2 - #1) + m(c p - c„)( d 3 - d 2 ) - mc v (d 4 - d 3 ) 



mc p (d 3 - 0 2 ) 



1 di/di — 1 
ft 0 2 ^3/^2 - 1 



This equation can be simplified by using the relationship for the isentropic process 
(Equation 2.70) and the relationship for the isobaric process (Equation 2.72). 
Additionally, the following relationship is used for the isochoric heat loss: 

d4 _ P4 _ P4P2 _ V 3 Vi 

di pi p 3 pi V 4 k V 2 k 

In that, we have p 2 = p 3. 



#4 _ fdA K/K_1 (d 2 \ K/K_1 = fdA K/K_1 \ K/K “ 1 

“ w UJ “UJ V0 3 ; 

This results in: d 4 /di = ( 03 / 0 2 ) K = p K , which yields the thermodynamic effi- 
ciency of the Diesel engine: 



Pth = 1 - 



1 1 p K - 1 

£ k_1 ft p — 1 



(2.73) 



It can be seen that the efficiency p t h decreases as the load p is increased. At high 
loads, the diesel engine has a lower efficiency compared to the SI engine, supposing 
the same compression ratio e for both (see Figure 2.13). The compression ratio 
for Diesel engines is however much higher than for SI engines to improve the 
thermodynamic efficiency. 



2.2.3 Seiliger Process 

The Seiliger process models the thermodynamic process in automotive engines 
much better than the previously described models of SI and Diesel engines. Fig- 
ure 2.12 shows that the combustion is now divided into two parts: In the first 
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V 3 > 



Vi 




Figure 2.12 pV-diagram of Seiliger process 



part, the gas is heated in an isochoric process between step 2 and step 3. In 
the second part the gas is expanded in an isobaric state change between step 3 
and step 3’. The cycle is characterized by the compression ratio e = V 1 /V 2 , 

the injection ratio p = Vy /V 3 , and the pressure ratio: 

X=~ (2-74) 

P2 

The different steps of the cycle are the following: 

1 — > 2 : Isentropic compression, dq = 0: 

dq = du + dw = 0 
91,2 = 0 

dw = —du=—mc v d , d 
w i j2 = -mc v (il 2 -i?i) 

2 — » 3 : Isochoric input of thermal energy, dV = 0: 

dw = p dV = 0 

tt>2,3 = 0 

dq = du = mc v dfl 
92,3 = mcyifi 3 -i? 2 ) 

3 — > 3' : Isobaric input of thermal energy, dp = 0: 

dq = dh — V dp = m c p dfl 
93,3' = m.c p (d 3 '-i? 3 ) 

dw = p dV = mRdD 
= mR(dz>- 1 ? 3 ) 



^ 3 , 3 ' 
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3' — * 4 : Isentropic expansion, dq = 0: 

93', 4 = 0 

dw = —du = —mc v d0 
W3',4 = -mc v (0 4 -1? 3 /) 

4 — > 1 : Isochoric heat loss, dF = 0: 

dw = p dV = 0 

W4 t i = 0 

dq = du + dw = mc v did 
(?4,i = mc v (0 ! - i? 4 ) 



The thermodynamic efficiency of the Seliger process is then: 



Vth 



U>1,2 + ^2,3 + W 33 ' + W3' )4 + W 4 ,i 
<72,3 



. -d 4 /-di - 1 

# 3 /t?i - fV'di + k(i? 3 '/^i “ ^3/^1) 
The isentropic process is characterized by: 



^2 _ K _i 

t?i 

which yields: 

— = 

0:i> \V 4 ) 

(VrViy - 1 = (v^vy 
\v :i vj \v 3 vj 

- (I)' -1 

In the isochoric process, the relationship 



$3 P3 
T = — = X 
02 P2 

can be used, and the temperature ratio in the isobaric process is given by: 

$3' _ V 3 > _ 

Therefore, the following temperature ratios may be expressed as: 



03 


03 02 


K — 1 


01 


02 01 


= X£ 


$3' 


03' 03 


K , — 1 


#1 


03 01 


= P X e 


04 


04 03 ' 


(P\ K ~ l 


01 


03' 01 


-\e) 
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Figure 2.13 Thermodynamic efficiency rjth depending on compression ratio e and load 
p (Equation 2.75) 



The thermodynamic efficiency of the Seliger process can be simplified to: 



_JL X p K - 1 

e K_1 X-1 + «X(P-1) 



(2.75) 



It can be seen that the thermodynamic efficiency of the SI engine is obtained 
when p = 1 and that of the Diesel engine when % = 1 (Figure 2.13). 



2.2.4 Comparison of Different Engine Concepts 

The in-cylinder pressure during combustion is plotted over the crankshaft angle in 
Figure 2.14. The compression ratio e for the SI engine is limited by the maximum 
allowable pressure ps during the combustion process. Under part load conditions, 
the maximum pressure of a cycle is far below this limit, since the SI engine power 
output is modulated by throttling the air intake thus modulating pi (Figure 2.15). 
A low compression ratio e is also helpful to reduce knocking of the engine and 
to meet material demands. In contrast to SI engines, the maximum pressure for 
Diesel engines is closer set to the maximum allowable pressure p 3 . As the Diesel 
engine is unthrottled (modulation of p) 1 it can afford higher compression ratios 
£ = V 1 /V 2 than the SI engine (Figure 2.16). 

The four-stroke engine works intermittently: A hot combustion cycle is al- 
ways followed by a cool gas exchange cycle. In the first cycle, peak temperatures 
of 2500 — 2800 °C occur. Compared to that, the maximum temperature of gas 
turbines must be kept much lower (1300 °C ), since the combustion process is con- 
tinuous. Diesel engines have a higher thermodynamic efficiency than SI engines 
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piston 




Figure 2.14 In-cylinder pressure over crankshaft angle during combustion process 



at low and medium power (Figure 2.13). At turbocharged Diesel engines also p\ 
is modulated. In order to stay below maximum temperatures, the compression 
ratio e is reduced at turbocharged Diesel engines. The absolute effective power 
output of combustion engines depends on the displacement volume V d . the spe- 
cific work output per power-stroke w e and the number of crankshaft revolutions 
n: 

TL 

P e = w e V d - (2.76) 

where: 

P e is the absolute effective power output 
w e is the specific work referred to the displacement volume V d 
V d is the displacement volume V\ — V 2 
(here, only one cylinder is regarded) 
n is the number of crankshaft revolutions per minute 

The factor 1/2 is necessary in determining the power output of four-stroke engines 
since only every second cycle contributes power. The mean piston velocity s 
depends on the number of crankshaft revolutions per minute n and the maximum 
piston stroke 2 r from top dead to bottom dead center. The mean piston velocity 
is: 



s = 4nr 



(2.77) 
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A typical value in car engines is s = 15 m/s at maximum power output. Equa- 
tion 2.76 can now be written as: 



Pe_ 

v d 




(2.78) 



Please note that the specific power P e /Vd is inversely proportional to the maxi- 
mum piston stroke 2 r (r is the crankshaft radius). Table 2.1 gives some examples 
of engine types and their characteristic values. It can be seen that the specific 
power decreases with increasing dimensions r at large engines. They are mainly 
used in ships because of their thermodynamic efficiency. When size and weight 
are major concerns at high power levels, gas turbines are preferred e.g. in air- 
craft. 



2.3 Alternative Combustion Engines 

2.3.1 Gas Turbine 



In the gas turbine engine fuel is combusted in a continuous process. The air has 
to be compressed before it flows into the combustion chamber. There the fuel 
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Figure 2.16 Load behavior of Diesel engine (Modulation of p), 1 bar = 10 5 Pa 



is injected and burned. The resulting expansion of the gas is used to turn the 
turbine shaft. Therefore, the effective power is the power of the turbine minus 
the power used for the compressor. The process can be modeled as a Brayton 
cycle (or Joule process) which is shown in Figure 2.17. The different steps of 
the Brayton cycle: 

1 — > 2 : Isentropic compression, dq = 0: 



dq = du + dw = 0 
Qi,2 = 0 

dw = — dh = — to c p did 

wi , 2 = -mcp(i? 2 - i?i) 



The kinetic energy w i 2 is provided by the compressor. 



2 — > 3 : Isobaric input of thermal energy, dp = 0: 



dq = dh — V dp = m c v di) 
<?2,3 = mcp(id 3 -i? 2 ) 

dw = p dV = mRd'd 

= ?ni?(d 3 -'d 2 ) 



^ 2,3 
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Table 2.1 Specific power output of various combustion engines 



Engine type 


s 

[m/s] 


w e 

[J / cm 3 } 


2 r 

[cm] 


Pe/V d 

[kW/dm 3 ] 


Car engine 


15 


1.6 


8 


75 


Truck engine 


11 


1.9 


14 


37 


Big four-stroke diesel 


9 


2.0 


32 


14 


Big two-stroke diesel 


7 


1.8 


60 


5 




Figure 2.17 pV-diagram (left) and i?S-diagram (right) of the Brayton cycle 



3 — > 4 : Isentropic expansion, dq = 0, generation of work: 



93,4 = 0 

dw = — dh = — m c p dd 

W3,4 = —mc p {d 4 -^ 3 ) 



4 — > 1 : Isobaric heat loss, 



dp = 0: 






dq = 


dh — V dp = to c p dd 


94,1 = 


mc p (di — 


d 4 ) 


dw = 


mRdd = 


0 


W4 , 1 = 


mR(di — 


d 4 ) 



The thermal energy < 3^1 and mechanical work u> 4 ,i are both negative be- 
cause of $1 < d 4 . 
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The thermal efficiency of the gas turbine is: 

92,3 + 94,1 

r\th = 

92,3 

mc p (il 3 — $2) + mc p (#i - $4) 
mc p {d 3 - 1I2) 

1 - #l/l?4 
$3 1 - ^2/^3 

As pointed out before, the isentropic process is characterized by: 

= my -1 1 _ 

$2 VvJ £K_1 ^3 

The relationship = ^2/^3 for the isentropic process and the injection ratio 

p = $3/ $2 are used to simplify the thermodynamic efficiency: 

nth = 1 - 1 {2.79) 
p v 2 

The efficiency of the gas turbine depends on the temperature ratio i?4/i?2- Rea- 
sonably high efficiency levels can only be reached, if the thermal energy of the 
outgoing air (74 4/ is used to heat up the incoming air in front of the compressor 
92,2'- Figure 2.18 shows the i9S-diagram supposing a complete recirculation of 
thermal heat. Assuming $2 = $4' and $2' = $4 , the necessary input of thermal 
energy for each cycle can be reduced to: 



92 ', 3 = mc p (il 3 - $ 2 ') = mc p (tl 3 - $ 4 ) 
and the heat loss is also reduced to: 

94', 1 = mc p {‘d 1 — $ 4 ') = mc p {d\ — r d 2 ) 

Equation 2.79 can be modified: 

92', 3 + 94', 1 

nth = 

92', 3 

$2 1 ^ ■& l/#2 
1 “ ^4/^3 

Finally, by using the same relationships as before, the thermodynamic efficiency 
of the gas turbine with complete heat recirculation is: 

mu = 1 - 1 (2-82) 

P 

The thermal efficiency depends only on the load p = i?3/i?2- The higher the load 
the better the efficiency. The maximum temperature $3 is mainly limited by the 
material from which the gas turbine is constructed. Equation 2.82 is not valid 
for extremely small loads as the heat recirculation does not work properly under 
these conditions. The main advantages of gas turbines are the following: 



(2.80) 

(2.81) 
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Figure 2.18 $S-diagram of the Joule process with heat recirculation 



• Gas turbines are much smaller and lighter compared to four-stroke piston 
engines: The weight per power output in gas turbines is 0.3 — 0.5 kg/kW 
and for four-stroke piston engines it is 1.0 — 1 .5kg /kW. 

• High torque can be generated even at very low revolutions. 

• Different types of fuel can be used for combustion: multi-fuel capability. 

• Gas turbines are easy to start even at low temperatures. 

• Low vibrations because of a continuous combustion and the rotary motion. 

• Long service intervals between required maintenance. 

• Reduced emissions of noxious exhaust gases. ECE test results: 

CO = 20 g, HC = 0.8 g, NO x =2 — 3 g per test. 

However, there are some disadvantages: 

• Low efficiency for low loads. 

• Poor dynamic behavior during transients. 

Gas turbines are mainly used in air planes because of their low weight. Their effi- 
ciency can be increased by raising the maximum allowable temperature. Ceramics 
like AI 2 O 3 or laminated silicon-carbon materials are used for the construction to 
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piston 



/ | \ 

heater regenerator cooler 



Figure 2.19 Mechanical representation of the Stirling engine 




allow higher temperatures. An example of temperature values for two different 
turbines types is given in Table 2.2. 

Table 2.2 Temperature values in °C at different locations of a metal and a ceramic 
gas turbine 



Location of measurement 


Metal gas turbine 


Ceramic gas turbine 


Compressor inlet 


230 


250 


Heat exchanger outlet (air) 


700 


950 


Combustion chamber outlet 


1000-1100 


1250-1350 


Heat exchanger inlet (gas) 


750 


1000 


Heat exchanger outlet (gas) 


270 


300 



2.3.2 Stirling Engine 

The Stirling engine is a piston engine which uses a continuous heat supply. It 
was invented by Robert Stirling in 1816. The cycle process has a high efficiency, 
comparable to that of the reference Carnot cycle. Even though the engine was 
redesigned by Philips Corporation (in 1938) and recently by other companies, 
some mechanical problems remain unsolved. Figure 2.19 shows the mechanical 
representation of the Stirling engine. It consists of three main parts: 

H heater the gas is heated from outside 

R regenerator the regenerator stores thermal energy 
C cooler the gas is cooled from outside 

The heater, regenerator and cooler are located in the middle. The pistons on the 
left and right side are linked mechanically. The cycle consists of the following 
four steps that can be seen in Figure 2.20: 
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1 — > 2 : Isothermal compression: $ i = 1)2 

The heat is absorbed by the cooler when the two pistons move to the left 
to compress the gas. 



w 1,2 

91,2 



mR'di In 



Vi 



Wl,2 



The emission of heat < 71,2 is equivalent to the input of kinetic energy u>i i2 - 

2 — > 3 : Isoclroric input of thermal energy: V 2 = V 3 

The two pistons move simultaneously to the left and the gas is heated by 
the regenerator. Both, pressure and temperature are increased. 



<72,3 = mc^d 3 - d 2 ) 



3 — ■> 4 : Isothermal expansion: $3 = $4 

The thermal energy is supplied by the heater, when the pistons move back 
to the right. 



W 3 ,4 

93,4 



V 4 

mRi93 In— 
V3 

W3,4 



The mechanical work w 3t4 is equivalent to the thermal energy q 3 4 . 
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Figure 2.21 The four steps of the Stirling cycle 



4 — » 1 : Isochoric heat regeneration: V 4 = V\ 

As the pistons move simultaneously to the right, thermal energy is stored 
in the regenerator. The pressure as well as the temperature drop to lower 
levels. 



94,1 = mc v ( - $ 4 ) 



By exploiting the fact that $1 = $2 and $3 = $4 for the isothermic processes, 
and V] = V 4 and V 2 = V 3 for the isochoric processes, the thermal efficiency of 
the Stirling engine is: 



Vth 



91,2 + 93,4 _ 1 $1 

93,4 $3 



(2.83) 



Hence, the efficiency depends only on the temperature ratio il i/Ds- For example, 
a temperature ratio of 80°C'/600°C' results in a thermodynamic efficiency of: 



Vth = 0.87 



(2.84) 



which is very close to the efficiency of the Carnot cycle. The main advantages of 
the Stirling engine are: 

• Engine is independent of the heat source. Instead of combusting fossil fuels, 
alternate heat sources such as solar heat could be employed. 

• High (theoretical) efficiency. 

• Very quiet. 

• Reduced emission of noxious exhaust gases. ECE test results: 

CO = 4 - 69 , HC = 0.5 — 2g, NO x = 0.6 - 2.0 g per test. 

O 11 the other hand, there are some disadvantages: 

• Expensive to construct. 

• Regenerator: conduction and storage of heat are difficult to combine. 
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Figure 2.22 Steam engine 



• Heat resistant materials needed. 

• A heat exchanger for the cooler is needed to increase the efficiency. It 
increases however volume and costs. 

Experimental Stirling engines with temperatures of 40 — 80 °C/600 — 650 °C can 
reach an effective thermodynamic efficiency of 

r} e ff = 0.35-0.40 (2.85) 

which is much lower than the theoretical value (Equation 2.84). 



2.3.3 Steam Engine 

The steam engine introduced by James Watt is the oldest engine using continuous 
combustion. The steam is transferred from the boiler to the cylinder. The piston 
is moved by the expanding steam. The linear movement of the piston is translated 
into a rotation of the crankshaft by the connecting rod. The different steps are 
illustrated in Figure 2.22 and Figure 2.23. 

The different steps of the steam engine cycle are: 

1 — > 2 : Isoclroric compression, followed by an isothermal expansion: 

This process can be divided into two parts: First, hot steam is injected into 
the cylinder through the open valve at constant volume (1 — > 1') where 
Vy = V\. Second, the gas expands at a constant temperature (1' — > 2) 
where fly = $ 2 : 



W 1 , 1 ' 
Ql,V 

Wl',2 

Ql>,2 



0 

mc v (fly — di) = mc v (fl 2 — i?i) 

V 2 V 2 

m R fly In — = m.Rfl 2 In — 
Vy V± 



wy, 2 
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Figure 2.23 pV-diagram of the steam engine 



2 — > 3 : Isentropic expansion: 

After the valve is closed, the expansion is continued until the maximum 
volume is reached. 

92,3 = 0 

W2,3 = mc v (tf 2 -# 3 ) 

3 — > 4 : Isochoric heat regeneration and isothermal compression: 

This process can be divided into two steps: The pressure drops at constant 
volume after the valve is opened at 3 — » 3' where V 3 = V 3 >. Second, the 
steam is compressed isothermally at 3' — > 4 where $ 3 / = $ 4 . 



W3, 3' 


= 0 




93,3' 


= — 1 ? 3 ) 


= mc v (il4 - i? 3 ) 


W3',4 


v 4 

= m R, dy In — 

V 3 , 


V 4 

= m.Rili In — 

V3 


93', 4 


= W3’A 





4 — » 1 : Isentropic compression: After the valve is closed, the gas is compressed 
mechanically: 

94,1 = 0 

w 4,1 = mc v ( $4 — i?i) 

The mechanical work 104,4 is negative. 
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The thermal efficiency of the steam engine is expressed by: 

Wl',2 + U>2 , 3 + w 3 ',4 + 104,1 

nth = — 

91,1' + 91', 2 

_ 1?3 - #4 + (k - 1)^4 ln(y 3 /V4) 
"^2 — di + (k — l)i?2 ln(V2/Vi) 



By inserting the compression ratio e = V 3 /V 1 , the partial compression ratio 
p = V 2 /V 1 = V 3 /V 4 and the pressure ratio % = pv/pi, this leads to: 



nth = 

An example is given for p = 



1 1 (k - 1)(1 + In p) 

e*" 1 (X- 1) + («- l)x ln d 
2 and \ = 10 : 



, 1.065 

’M = 

which is rjth = 0.31 for a k = 1.4 and e = 3. 



( 2 . 86 ) 



(2.87) 



Advantages of the steam engine: 

• Engine is independent of the heat source: multi-fuel capability. 

• Noxious exhaust emissions are low because of continuous combustion. 

• High torque at low revolutions. 

Disadvantages of the steam engine: 

• Heavy weight. 

• Poor thermodynamic efficiency. 

• The water in the boiler needs to be heated before the engine can be started. 



2.3.4 Potential of Different Fuels and Propulsion Systems 

Table 2.3 illustrates that a constant amount of stored energy varies considerably 
in its volume and weight. Standard lead batteries are much too heavy. Other 
types of batteries are lighter, but they are still not comparable to the weight of 
ordinary fuel. Power is dissipated in the charging and discharging processes of 
the battery, reducing the overall efficiency. Eventually, battery driven vehicles 
with a reduced buffer size may be used in special applications at short distances. 
Another promising approach is that of hybrid vehicles, where an internal com- 
bustion engine is combined with an electrical motor. The electrical motor may 
be activated to smooth out transients of the combustion engine and the drive- 
line, contributing to reduced noxious emissions. Under partial load conditions 
the combustion engine can also charge the battery, so that battery volume and 
weight are significantly reduced. Hydrogen H 2 gas is too voluminous to be used 
as adequate energy source. It can be stored either at an extremely cold tempera- 
ture of 20 K or at relatively high pressure at room temperature. Over long time 
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Table 2.3 Typical storage volumes and weights of different energy sources with an 
energy of 1000 kWh. 



Source 


Volume 

V in [Z] 


Mass 
mi in [kg] 


Tank 

m 2 in [kg] 


Mass + Tank 

mi + m 2 in [kg] 


Fuel 


117 


83 


21 


104 


Diesel 


102 


85 


17 


102 


Methanol 


224 


180 


41 


221 


Liquid gas 


153 


78 


90 


168 


Methane 


259 


72 


500 


570 


H 2 , liquid 


426 


30 


142 


172 


H 2 , hydride buffer 


200 


30 


970 


1000 


Battery (lead) 


5000 


0 


10000 


10000 



energy demand 

[kWh/ 100 km] H fuel production 

A H driving operation 
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Figure 2.24 Energy demand of different engine concepts [74] 



periods, H 2 leaks through even thick walled steel tanks. In hydride buffers, H 2 
is chemically bound. Since hydrogen burns at high combustion temperatures, 
emissions of nitrogen oxide ( NO x ) become a problem. 

Fuel cells produce electrical energy directly at low temperatures. Thermal 
efficiencies of 70% are reached for the synthesis of H 2 and 0 2 . The storage of 
hydrogen is again the problem. If H 2 must be however generated from natural 
gas or from methanol, efficiencies become much lower. The task is to generate 
the exact amount of hydrogen from e.g. methanol even under realtime transient 
engine conditions. For this the fuel conversion process can be modeled, and the 
actual masses reacting in the conversion process be estimated in realtime, as a 
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basis for state space control. Fuel cells appear to be a promising alternative to 
combustion engines. In Figure 2.24, the relative energy requirements to move a 
vehicle by 100 km are shown for different propulsion systems. 




3 Engine Management Systems 



3.1 Basic Engine Operation 

3.1.1 Effective Work 



Four-stroke engines are characterized by two alternate cycles: In the first cycle, 
equivalent to the first and second strokes, the gas is compressed, combusted and 
expanded. In the second cycle, equivalent to the third and fourth strokes, the gas 
is transferred to the exhaust pipe and the cylinder is filled with fresh air from the 
intake manifold. Figure 3.1 shows the two cycles. The crankshaft is turned 360° 
per cycle. SI and diesel engines are controlled differently: In diesel engines, fuel 
is directly injected into the combustion chamber. The amount of injected fuel 
per stroke is then proportional to engine torque. The amount of air is almost 
constant at a given speed. In SI engines, the amount of fuel as well as air is 
controlled. When the fuel is injected into the intake manifold, a homogeneous 
air-fuel mixture is sucked into the cylinders. The mechanical work generated in 
the combustion cycle can be obtained by integration in the pE-diagram. The 
mechanical work can be normalized when relating it to the displacement volume 
V d : 




where: 



Vd = CYL ■ (Ei — V 2 ) is the displacement volume of all cylinders 
CYL is the number of cylinders 

Wi is the (normalized) indicated specific work. 
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Figure 3.1 pV-diagram of four-stroke combustion engine 



The value of uj, can be determined by measuring the in-cylinder pressure during 
a cycle. An indicated specific work of 1 J/cm 3 is equivalent to a mean pressure 
of p = 10 bar (= 10 6 Pa). The transfer of the combustion torque to the engine 
torque available at the crankshaft can be calculated from the following motion 
equations. 



The piston stroke from Top Dead Center (TDC) is 

s(acs) = i(l - cos/3) + r(l - cos a C s) 
From Figure 3.2 we get 



1 - , 2 sin a C s 



/sin/3 = r sin acs 
cos /3 = 

which yields the piston stroke as 

s(«cs) = r( 1 - cosacs + y(l-\/l-^ sin2 a ° s 



(3.2) 



(3.3) 



At Top Dead Center, we have acs — 0, s(acs) = 0, and at Bottom Dead Center 
acs = 7r> s(acs ) = 2r respectively. The derivatives of the piston stroke are 



ds 



sm a C s 



sin acs cos acs 



dacs 



r 



1 7*2 • 2 

1 — ~[t sm acs 



(3.4) 
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TDC 




s : piston stroke 

r : crankshaft radius 

cues '■ crankshaft angle 

l : connecting rod length 



Figure 3.2 Piston and crankshaft motion 



and 



d 2 s 

da 2 cs 



cos acs 



fcos 2 acs - sin 2 acs ) + tj sin 4 acs 



\ 



V 



• 2 

- w sin acs 



. (3.5) 



7 



These derivatives with respect to crankshaft angle can be related to the deriva- 
tives with respect to time as follows: 



s = 



ds 


ds 


dacs 


ds 








dt 


dacs 


dt 


dacs 


acs 






d 2 s 


d f 


ds 


da C s\ 


d f 


ds \ 


dacs 


dt. 2 


dt \ 


da C s 


dt ) 


dt \ 


da C s J 


dt 



dacs dt 2 



d 2 s 
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The indicated specific work can be written as 

. CYL 

l 

V d 



Wj = 



Y & '^2(pj(a C s) - Po) A p dS ^ CS ^ dacs 

3 = 1 

— T comb {acs)dacs ■ 



V d 

The combustion torque at the crankshaft is thus defined as 

CYL 



T CO mb(acs) = ^2 ( Pj( a CS ) - Po) A p 



3 = 1 



dsj 

dacs 



The piston strokes in different cylinders are shifted by phase. 

47T 



Sj(acs) = s a C s - (j ~ 1) 



CYL 



j = l,...,CYL 



(3.7) 



(3.8) 



(3.9) 
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rrif ■ Hf combustion heat 

Qhi,th theoretical heat, loss 
Wth theoretical work 

heat loss caused by 
yhi,r incomplete combustion 

Wi,hp high pressure work 

y vJi,ip low pressure work 

Wi indicated work 

w f r frictional work 
w e effective work 

Figure 3.3 The effective work delivered by the engine is much lower than the thermal 
energy caused by combustion. 






The average combustion torque is 



I'rornh 




Pi 

dies 



T CO mb{acs)dacs 



(3.10) 



where Pi is the mean indicated power. The total indicated work w, : V d can now 
be written at stationary engine operation as 



WiV d = 47rT comfc = 47 t- 



Pi 

&CS 



MrP 

2irn 



2P 

n 



and the normalized work 

= • < 311 > 

where n = acs/{ 2tt) is the engine speed. In reality, the effective work w e per 
volume is much lower than the indicated work Wi (see Figure 3.3). The effective 
thermodynamic efficiency r) e is at constant fuel flow 



Pe _ W e V d n 2 _ W e V d 

rh f H f ~ 2m f nH f ' CYL ~ m f H f ' CYL 



(3.12) 



where: 

P e is the effective power in W 
w e is the effective specific work per cycle in J/cm 3 
mf is the mass of fuel measured per cylinder in kg 
rhf is the fuel flow in kg/s 

Hf is the specific energy of the fuel released in the combustion J/kg 
V d is the total displacement volume in dm 3 

{ Vd/CYL displacement volume per cylinder) 




3.1. BASIC ENGINE OPERATION 



47 



The indicated thermodynamic efficiency (friction not considered) is: 



Wj V d 
2 m f H f ' CYL 



(3.13) 



Some examples of typical values for the indicated efficiency are given in Table 3.1. 



Table 3.1 Indicated specific work Wi, theoretical heat loss qhi,th, and realistic heat 
loss qhi, r for different engine types, related to fuel combustion heat. 



Engine Type 


SI 


Diesel 


Big Diesel 


Wi 


33-35 % 


40-43% 


45-48 % 


Qhl,th 


23-28 % 


22-25 % 


12-14% 


Qhl,r 


37-44% 


35-40 % 


26-33 % 



3.1.2 Air- Fuel Ratio 



The ratio of air to fuel is very important for the combustion process of inter- 
nal combustion engines. There are several effects that have an impact on the 
amount of air m a transferred to the cylinder: Throttling of the air flow by the 
throttle butterfly, aerodynamic resistance and resonances in the intake manifold, 
rebounding of already burned gases from the cylinder into the inlet pipes and 
other effects. The amount of air which would theoretically fit into a displacement 
volume Vd under the normalized pressure po = 1.013 bar and the normalized air 
density po = 1-29 kg/m 3 is expressed by m a ^h = PoVd- The ratio of real to 
theoretical value is equivalent to the relative air supply: 



m a 

TTl a,th 



(3.14) 



Similarly, the ratio of measured fuel mass m/ to theoretical fuel mass rtif^h is 
equivalent to the relative fuel supply: 



X f 



m f,th 



(3.15) 



The theoretical fuel mass mf^h is equivalent to the mass of fuel needed for an 
ideal stoichiometric combustion with the oxygen. Under normal conditions the 
stoichiometric ratio for gasoline is: 

L st = nl,hth = 14.66 (3.16) 

m f,th 

The air-fuel ratio lambda is defined as: 
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It can be extended: 



m a inf,th _ 1 m a 
mf m a ,th L s t rn f 



(3.18) 



For an ideal stoichiometric combustion, this ratio is equivalent to one: A = 
1. The air- fuel ratio has an impact on the effective work w e and the effective 
thermodynamic efficiency r] e . The air- fuel ratio can be influenced in two different 
ways, by variation of A a or of A/: 



1. Variation of A/ at a given A 0 : 

Typical applications are SI engines operating around a stoichiometric air- 
fuel ratio. The relative air supply A a is determined by the driver. 



Lean operation (A > 1): Less fuel is injected than needed for stoichio- 
metric combustion (reduced A/). Due to a reduced high pressure work 
Wi^hp the effective work Wi t h P decreases. In the range of 1 < A < 1.1, 
the thermodynamic efficiency rj e increases however, caused by higher 
combustion peak temperatures. This results in high emissions of ni- 
trogen oxides NO x . If A is further increased, rj e will decrease because 
of an even lower high pressure work wi t h P at a given low pressure work 

W i,lp • 

Rich operation (A < 1): More fuel is injected than needed for stoichio- 
metric combustion (higher A/). The fuel surplus increases both high 
pressure work Wij l/p and effective work w e . Below A < 0.9 incom- 
plete combustion results in high emissions of hydrocarbon HC in the 
exhaust gases and in a decreasing effective work w e . At A < 1 the 
thermodynamic efficiency r) e is always decreased. 



2. Variation of A a at a given A/: 

Typical applications are lean-burn Si-engines at part load and Diesel en- 
gines. The relative fuel supply A f is determined by the driver. 

Lean operation (A > 1): Alore air is admitted than needed for stoichio- 
metric combustion (increased A a ). Therefore, high pressure work Wi t h P 
is increased while low pressure work Wij p remains constant. Both effec- 
tive work w e and thermodynamic efficiency r] e are increased. It should 
be mentioned that lean gas mixtures are less flammable. In Si-engines 
delays between spark ignition and complete combustion increase. 
Misfiring at SI engines must be avoided by e.g. direct injection of 
a fuel stratified charge into the cylinder, which forms an enriched 
mixture around the spark plug. This operation is quite similar to 
that of Diesel engines. Combustion is either triggered by a spark or 
by self-inflammation due to high compression ratios. The engine can 
only operate up to its maximum gas load (maximum A a ). Due to lean 
operation its maximum power according to the displacement volume 
is, however, not reached. 

Rich operation (A < 1): Less air is admitted than needed for stoichio- 
metric combustion (decreased A 0 ). This leads to a decrease in both 
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(1) SI engines at maximum, power output 

(2) Stochiometric SI engines 

(3) Diesel and lean — burn engines 

(4) Moderately lean SI engines 



de , [%] 




Figure 3.4 Effective work w e and thermodynamic efficiency r] e of combustion engines 
depending on variation of A a or A/. 



efficiency r] e and effective work w e . Incomplete combustion results in 
higher hydrocarbon HC emissions and in a reduction in high pressure 
work Wi t h p - 

Figure 3.4 shows the dependency of the effective work and effective thermody- 
namic efficiency over A, assuming an optimal control of fuel injection and ignition 
timing. 

Engines may use recirculated exhaust gas instead of fresh air to increase the 
relative air supply A 0 . As long as sufficient air is available for the combustion 
this is similar to a higher air- fuel ratio A. Exhaust gas recirculation reduces the 
emission of NO x due to lower combustion peak temperatures. 

3.1.3 Engine Concepts 

The SI engine is controlled by the relative air supply A a . This is done by throt- 
tling the air flow into the engine. The relative fuel supply A/ is subsequently 
regulated by the engine management systems to maintain the desired air-fuel 
ratio A. The range of A is limited by the ability to inflame air-fuel mixtures by 
spark ignition. Conventional SI engines operate on approximately homogeneous 
mixtures (0.9 < A < 1.3). Lean-burn engines operate at very lean mixtures equiv- 
alent to Diesel engines. Combustion is ensured by directly injecting a stratified 
charge of rich air-fuel ratio around the spark plug. 

The Diesel engine is controlled by the relative fuel supply A/. The intake 
manifold is not throttled. The relative air supply A a is always at its maximum. 
Therefore, the air-fuel ratio A changes within a large range. The inflammation 
of extremely lean mixtures is still possible because of the non-homogeneous fuel 
distribution in the combustion chamber. Such inhomogeneous mixtures burn with 
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w e w e 




Figure 3.5 Effective work w e over engine speed for SI engines (left figure) depending 
on the throttle angle at and Diesel engines (right figure) depending on relative fuel 
supply A f . 



a yellow flame. The average air- fuel ratio should not be below A = 1.3 to avoid 
the generation of too much soot. Since the effective work w e is basically given by 
the amount of injected fuel, the fuel supply must be cut-off when reaching the 
maximum engine speed. Otherwise the engine power would continue to increase 
with speed resulting in a self-destruction of the engine. Fuel may be injected 
in several steps. A first small amount of fuel starts the combustion process 
more smoothly. The second main injection then results in lower peak pressures 
and temperatures, yielding lower NO x emissions and less combustion noise. An 
injection at the end of the combustion cycle heats up the exhaust pipe and 
exhaust gas treatment systems such as a soot filter. When Diesel engines are 
operated with very high exhaust gas recirculation rates, then they generate few 
noxious raw emissions [28]. 

Lean-burn SI engines are a compromise between diesel and stoichiometric 
SI engines. Driving at part load, the air-fuel ratio is very lean. By properly de- 
signing injection pressure, spray cone and air turbulence, an enriched stratified 
charge is assembled around the spark plug. The resulting combustion is equiv- 
alent to that of Diesel engines (inhomogeneous mixture, yellow flame). Lower 
noxious emissions can be achieved by separating into e.g. two injections. First, 
about 1/4 of all fuel is injected, which forms a lean homogeneous mixture in the 
combustion chamber. This lean mixture is also less sensitive to knocking. Sec- 
ond, 3/4 of fuel is injected as a stratified charge. After burning this rich mixture, 
the homogeneous lean mixture burns in a blue flame, resulting in a reduction of 
noxious emissions. At high loads, the operating conditions are shifted from very 
lean to stoichiometric mixture. At a given engine displacement volume, more 
fuel can be combusted from stoichiometric mixtures, increasing power output. 

Figure 3.5 shows the dependency of the effective work w e over engine speed for 
SI and Diesel engines. The displacement volume of a naturally aspirated diesel or 
lean-burn SI engine must be 60 % higher than that of a stoichiometric SI engine 
to obtain the same maximum power output. Therefore diesel and lean-burn SI 
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Vi $ isentropic 




Figure 3.6 Isentropic temperature Disentropic and pressure p over time. The conditions 
for self inflammation (a, b, c ,d) are explained in the text below. 



engines are often turbocharged which increases the relative air supply A a at a 
given displacement volume Vd- 

3.1.4 Inflammation of Air-Fuel Mixtures 

The kinetic gas theory describes gases as a cloud of molecules with a given velocity 
distribution according to their temperature. The collision of different molecules 
will start a chemical chain reaction if their kinetic energy is over a certain ac- 
tivation energy E. The relative amount of effective collisions A is expressed by 
the Arrhenius law 

A = e ~ E/m . (3.19) 

The activation energy E is low for radicals (not saturated molecules). The prob- 
ability for a collision is increased by the concentration of molecules and by the 
temperature. A chemical reaction must be started by a high temperature. During 
the chain reaction, more radicals are generated than destroyed. Under appropri- 
ate conditions, a spontaneous spark ignition is sufficient to start the combustion 
process at the location of the spark plug. The air-fuel mixture must be within a 
certain range (0.9 < A < 1.3) and its pressure (or respective temperature) must 
be over a threshold for a certain period of time. The gas is compressed isentropi- 
cally, neglecting heat conduction. Figure 3.6 shows how pressure and temperature 
courses over time influence self inflammation under different conditions. 

1. In curve a, pressure and temperature rise immediately. In curve 6, inflam- 
mation starts only when reaching a higher pressure level than in curve (a). 
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Tid, 




Figure 3.7 Delay of self inflammation T,d over temperature $ for different pressure 
values p. 



In the case of (b), the increase in pressure and temperature was first delayed 
until tg. 

2. In curve c, temperature and pressure rise immediately, but level off at a 
lower level. Inflammation occurs after a longer time delay. In the case of 
(d) the rise is first delayed until time to- Self inflammation happens later 
than in (c). 

It can be seen, that self inflammation depends on something like the integral of 
pressure or temperature over time. Many parameters have an impact on the time 
delay of self inflammation, like location of the spark plug within the combustion 
chamber, etc. Woschni [130] gives an empirical formula of the inflammation delay 
time, depending only on mean temperature i9 and mean pressure p , without any 
integral portion. 

T id = 0.44 e 4650K /& (ILj (3.20) 

The inflammation delay T t( j over temperature i9 with the pressure p as parameter 
can be seen in Figure 3.7. The time delay r,; c | is reduced for high temperatures 
and high pressures. This is why e.g. turbo-charged Diesel engines operate with a 
start-of-injection angle which is approximately 10° (crankshaft angle) later than 
that of naturally aspirated Diesel engines. 

3.1.5 Flame Propagation 

The velocity Vfj of flame propagation depends on two components: 

1. Combustion velocity: The combustion propagates through the gas mix- 
ture, for example with a velocity of 1 m/s. 

2. Transport velocity: The burning gas itself is swirled as the rising piston 

generates turbulence in the combustion chamber. The transport velocity 
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Vf U [m/s] 




v fh [m/s] 

^ (a) static gas at po 

( b ) static gas at p > po 
2 - (c) turbulent gas at p > 




Figure 3.8 Flame propagation velocity Vfi over average piston velocity s (left) and 
over air- fuel ratio A (right). 



is approximately proportional to the piston velocity s which depends on 
the engine speed n. At low engine speeds, the transport velocity can be 
increased by a swirl inlet port generating a turbulent gas flow. The turbu- 
lences accelerate the combustion speed proportional to engine speed. 

Figure 3.8 shows the flame propagation velocity Vfi depending on piston veloc- 
ity s (left) and air- fuel ratio A (right). The inflammation delay time must be 
considered in the engine control to position the combustion process right over the 
downward moving piston. The time delay Tid must be convoluted to a crankshaft 
angle delay, increasing with the engine speed. The ignition angle on must there- 
fore be advanced over engine speed. 

Contrary to that, the position of the combustion process over the crankshaft 
angle is almost constant. This is due to the fact, that the flame propagation 
velocity at combustion is mostly determined by the transport velocity. Thus 
the angle position of the combustion process is independent of the engine speed. 
With higher engine speeds, flame transport velocity is increasing, speeding up 
combustion over time, leaving it however constant over engine speed. 

3.1.6 Energy Conversion 

The in-cylinder pressure can be plotted over time or over crankshaft angle acs- 
The angle of 360° is equivalent to a complete high pressure cycle. Figure 3.9 
shows the in-cylinder pressure over crankshaft angle. 

The gas is compressed by the piston in an approximately isentropic process. 
With ignition at on, the pressure rises only after time delay Tid- The maximum 
pressure varies from cycle to cycle. The inflammation delay Tid depends on tem- 
perature, pressure, air-fuel ratio and self inflammation time as described in the 
previous section. It also depends on the type of fuel being used. Figure 3.10 
shows some inflammation delays for different fuels over temperature. Oil compa- 
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P, [ bar ] 




Figure 3.9 In-cylinder pressure p over crankshaft angle acs- 



nies adapt their fuel to weather conditions (summer, winter). 

Turbulence caused by the upward moving piston has no impact on the time 
delay T lr i- For a correct ignition angle, this delay must be considered. The 
time delay is convoluted to an angle delay, increasing proportional to engine 
speed. Contrary to that, the engine speed has almost no impact on the position 
of energy conversion as turbulences increase the transport velocity with higher 
engine speeds. 

The energy conversion caused by combustion is shown in Figure 3.11 for 
different air-fuel ratios A. In these curves, the isentropic pressure curves are 
suppressed. The differential output of thermal energy per angle dQ/dacs (its 
gradient) is normalized to the total combustion energy Q com b ■ The shape of the 
relative energy conversion is therefore almost constant. 

If the air-fuel ratio is increased e.g. to A = 1.2 as shown in Figure 3.11, the 
ignition delay will rise. At a constant inflammation angle an the energy 
conversion is then retarded. Therefore, the ignition angle must be advanced to 
ai 2 , to compensate for the increased delay. The energy conversion returns to its 
previous position. It should be mentioned that a high air-fuel ratio A increases 
also the variance of the time delay r^. 

The ignition angle a,; depends on A which can be seen in Figure 3.12. The 
angle is computed by averaging the energy conversion over 0.1 %, 1 %, 10 %, 50 %, 
90 % points. The angles for «qi% and higher are almost independent of the air- 
fuel ratio A. In-cylinder pressure measurements can be used to control the ignition 
angle in a closed loop to maintain a constant position of energy conversion as 
shown in Figure 3.13. The angle of maximum pressure gradient max(dp/dacs) 
may be used as a control variable. The controller time constant must be relatively 
large because of the high delay time variances between consecutive cycles. Thus 
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Figure 3.10 Inflammation delay over temperature for different fuels. 
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Figure 3.11 Normalized energy conversion caused by combustion for different air- fuel 
ratios A. 
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Figure 3.12 Angle acs of energy conversion over air-fuel ratio A during ignition (left) 
and combustion (right) process. 



closed loop ignition control may be too slow for the dynamic response of the 
engine. 

The ignition angle is determined to find a compromise between fuel consump- 
tion, emissions or knocking (see Section 3.2.8). An equivalent procedure can be 
found for the fuel injection angle at Diesel engines. 



3.2 Engine Control 

3.2.1 Emissions of Internal Combustion Engines 

Mixture formation can be achieved by manifold or by in-cylinder injection. With 
sufficient time the mixture is distributed homogeneously in the cylinder with an 
air-fuel ratio in the range of 0.9 < A < 1.3. For very lean mixtures A > 1.3, a rich 
stratified charge must be concentrated in a portion of the combustion chamber. 

The combustion process is started by an electric spark at SI engines and by 
self-inflammation at Diesel engines. The inflammation is delayed as described in 
the previous section. 

• Homogeneous mixture, stoichiometric air-fuel ratio: The flame has a char- 
acteristic blue color. Almost no soot (carbon particulates) is produced. 

• Stratified charge, lean air-fuel ratio: The flame has a characteristic yellow 
color. Soot is produced. 

• Inflammation starts combustion from one location. 

The inflammation process depends on pressure p, temperature ■$, air-fuel ratio A 
and activation energy E of the fuel. For A < 1 the exhaust gases are generated 
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Figure 3.13 Closed-loop control of ignition angle at to maintain a constant position 
of energy conversion. 



according to the concentration ratio 

k= n C o-n H2 o (3.21.) 

nco 2 ’ h h 2 

This ratio is temperature dependant. A typical value for •& = 1850 °K is k = 3.6. 

The pollutant emissions like CO, HC, NO x depend strongly on the air-fuel 

ratio which is shown in Figure 3.14. 

A < 1: Increased emission of hydrocarbon HC and carbon monoxide CO. 

A = 1: Stoichiometric combustion. Very low emissions after three way catalytic 
converter. 

A « 1.1: Highest nitrogen oxide NO x emissions due to highest combustion peak 
temperatures. 

A > 1.1: Decreasing nitrogen oxide NO x concentration and lower combustion 
temperatures. Increasing hydrocarbon HC emissions at eventual misfires. 

A > 1.5: Lean operation. For very low emissions, a NO x reducing catalytic con- 
verter is required. 

The concentration of oxygen O 2 in the exhaust gas can be used to determine the 

air-fuel ratio A for A > 1 using a lambda-sensor. 
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Figure 3.14 Raw emissions of CO , HC, NO x and O 2 over air-fuel ratio A for SI 
engines. 

3.2.2 Fuel Measurement 

The air-fuel ratio A is an important variable for fuel control which is based on 
different control concepts: 

rich mixture A < 1: Maximum power output per displacement volume because 
of increased relative fuel supply Ay. It was used at high engine loads until 
1970. Nowadays it is only used for cold engines during the warm-up phase. 
High noxious emission rates. 

stoichiometric mixture A = 1: Acceptable power output. This ratio is re- 
quired for proper operation of three-way catalytic converters. At high en- 
gine loads, a good compromise between power output and exhaust emissions 
is achieved. 

moderately lean mixture 1 < A < 1.5: Good efficiency because of increased 
air supply A a , but high emissions of NO x . This method was used at part 
loads until 1980. 

lean mixture A > 1.5: High efficiency because of high A a . NO x emissions are 
still high, so that catalytic converters for NO x reduction are required. This 
method is used in lean-burn engines at part loads and in Diesel engines. 
Maximum engine power output cannot be reached. 

The reference torque desired by the driver controls either the relative air supply 
A a via the throttle angle a t at SI engines or the relative fuel supply A/ at Diesel 
engines. The amount of fuel being mixed with the air is regulated by the fuel 
control system to obtain a predefined air- fuel ratio A. There are two different 
injection systems: 
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1. Manifold injection: The fuel is injected into the individual inlet pipes in 
front of the inlet ports. There is at least one inlet valve for each cylinder. 
Problems may occur at idling because of incomplete fuel evaporation due to 
a low air flow velocity into the cylinders. Additionally, the distribution of 
air flow into the different inlet pipes may vary. The amount of injected fuel 
is less accurate at idling because electromagnetic injection valves are time 
controlled: Errors due to different rise and fall times have a larger impact 
on the amount of fuel injected at small injection times. The advantage of 
manifold injection is the creation of a homogeneous fuel distribution in the 
cylinder at A = 1. This burns with a blue flame. There are few restrictions 
for the design of the inlet pipes. The exchange of gas can be optimized 
without major effects on the injection system as it is located at the inlet 
valves. The inlet pipes are designed to produce acoustic resonances at 
low engine speed. This increases the relative air supply A a and the effective 
work w e already without turbo charging. The inlet valves are cooled by the 
evaporating fuel. The reduced gas temperature lowers knocking and allows 
higher compression ratios to increase efficiency. The injection timing is 
phase-shifted for each cylinder. It is aimed to terminate the injection just 
before the inlet valve is opened to avoid the emission of soot. The injection 
can be controlled individually and fuel supply can be individually cut off 
for each cylinder: limitation of engine speed and vehicle speed, fuel cut-off 
at coasting or cylinder switch-off at multi-cylinder engines. 

2. In-cylinder injection: The fuel is injected directly into the cylinders. At 
stoichiometric operation, a higher compression ration can be realized e.g. 
12 . . . 13, due to the cooling effect of the evaporating fuel. At lean operation, 
the aim is to assemble a sufficiently rich mixture (e.g. stoichiometric) in 
a limited portion of the combustion chamber, e.g. at SI engines around 
the spark plug at the time of ignition. The amount of fuel, the injection 
pressure (thus fuel atomization), the injector spray angle (width and depth 
of injection) and the injection timing are adjusted in each engine operation 
point. The swirl is controlled by special geometry of the piston head. 

At part load, the aim is to burn very lean air-fuel mixtures. This can be 
done by multiple injection. By an early injection, a homogeneous lean mix- 
ture is developing throughout the combustion chamber until the time of 
ignition, due to the swirl. The follow-up injection creates the richer strat- 
ified charge, which is burning fast with reduced cycle-to-cycle variations. 
Due to the reduced amount of the fuel in the stratified charge, peak tem- 
peratures are also reduced. The early generated lean homogeneous mixture 
is burned afterwards. Since the mixture is very lean, it takes significant 
time, limiting maximum allowable engine speeds. In the combustion of the 
stratified charge, soot was generated. In the blue flame of the homogeneous 
mixture, this can be effectively reduced. 

The total amount of injected fuel depends on the following parameters: 

• aspirated air flow per time m a 

• intake manifold pressure p m at SI engines 
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• throttle angle at and its derivative at SI engines 

• engine speed n 

• crankshaft angle acs and TDC signal of a reference cylinder 

• engine temperature d e 

• ambient air temperature d a 

• battery voltage Ub (indirectly) 

Major functions of the fuel control: 

• Control of injected fuel per time rhf, following the aspirated air per time 
m a , depending on the desired air-fuel ratio A. 

• Enriched fuel injection in warm-up phase of the engine after cold start at 
SI engines. 

• Increased relative air supply A a or relative fuel supply A / for the cold engine 
because of higher friction. 

• Compensation of intake manifold dynamics at SI engines. 

• Compensation of fuel film dynamics at manifold injection. This phenom- 
enon is also temperature dependant. 

• Fuel cutoff at coasting. This reduces overall fuel consumption by around 
5 %. 

• The measured air flow is eventually corrected for ambient air temperature 
d a and barometric pressure po changes. 

• Engine idle speed control. 

• Maximum engine speed limitation by fuel cutoff. 

• Lambda control of the air-fuel ratio. 

• Exhaust gas recirculation control. 

• Approximation of an engine torque signal, e.g. by means of a map for the 
thermodynamic efficiency r] e . 



3.2.3 Intermittent Fuel Injection 

Intermittent fuel injection has turned out to be more economical than continuous 
fuel injection, due to the different accuracies required in those systems. The 
power output of an engine varies within a wide range. Between idling power 
Pmin and maximum power P m ax > the ratio is about 1 : 100. 



Prr 



Prr 



= 100 



( 3 . 22 ) 
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The engine speed varies in the range of: 

“"max 



min 



10 



(3.23) 



At stationary engine operation, the amount of injected fuel per time to/ is pro- 
portional to the effective power output P e of the engine which is expressed in 
Equation 3.12. In this consideration air pulsation in the intake pipes are ne- 
glected. 

Supposing a continuous fuel supply to/, the relative error of the open loop 
fuel measurement at low loads should be 



A m 



<3% 



This will cause an absolute error related to maximum fuel flow 
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(3.24) 



(3.25) 



Hence, the absolute error at idling is 100 times smaller than the relative error. 
A continuous fuel injection system should be designed for the absolute error. It 
must be extremely accurate. High production costs are the consequence. 

Fuel measurement can alternatively be achieved by intermittent fuel in- 
jection. For each combustion cycle, a certain amount of fuel ?n/ is injected. 
The number of injections per second are proportional to the engine speed n. The 
amount of injected fuel per cylinder and combustion cycle is 



TO/ 




(3.26) 



where CYL is the number of cylinders of the engine. The factor 2 is due to 
the fact that air is combusted only every second cycle in the four-stroke process. 
Supposing a constant fuel flow to/ at stationary operation, the integration leads 
to: 



m > = it CYL . <“ 7) 

The ratio of maximum to minimum amount of injected fuel per cycle is given by 
Equation 3.12: 



= 10 

TFlmin ^ min ^ max 

Supposing again a relative error at minimum load of 



Ei 

Prr 



max 1 h min 



(3.28) 



Am 

min 



<3% 



(3.29) 



the absolute error related to maximum fuel per cycle is now 



Am 

TFlmax 



< 3 • 10~ 3 



(3.30) 



i.e. 10 times larger. Compared to a continuous fuel supply, intermittent fuel 
injection systems can be implemented with significantly lower accuracy require- 
ments and therefore can be produced at lower costs. 
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3.2.4 Injection Time Calculation 

The fuel supply is controlled by the injection time ti n j during which the injector 
valve is open. Therefore, the amount of fuel per injection into a cylinder can be 
calculated using the following relationship for constant air flow m a = const. 



_ m a _ 1 m a 2 

m/ “ L7t A “ L7t A V CYX 



(3.31) 



where L st = 14.66 (see Equation 3.16). The amount of injected fuel to/ is pro- 
portional to the injection time t,; n? and the square root of the pressure difference 
A p between fuel rail and manifold at manifold injection or between fuel rail and 
combustion chamber at direct in-cylinder injection [69]. The fuel density Qf and 
the effective opening area of the valve A e ff are assumed to be constant. 



TO/ 



Qf ■ A eSf ■ 




(3.32) 



At manifold injection, the pressure difference A p is around 5 bar. At direct in- 
cylinder injection, the pressure difference A p is up to 400 bar for SI and up to 
2000 bar for Diesel engines. 

The injection time at stationary engine operation is proportional to 



+ 1 rh a 2 

inj ~ A n CYL 



(3.33) 



and for a reference air-fuel ratio Ao a reference injection time to is proportional 
to 






1 ih a 2 
Ao n CYL 



(3.34) 



For arbitrary air-fuel ratios A we get 



tinj ~ 



(3.35) 



The injection time t ln j per combustion cycle depends on the following values: 

Mass air flow rh a : Must be measured. Systematic measurement errors at some 
sensors can be reduced by taking air density and temperature into account. 



Air mass per stroke ?n 0 : is computed by implementation of Equation 3.39. 

Reference air-fuel ratio Ao: Must be determined, e.g. stoichiometric. A look- 
up table can be implemented to compensate for possible errors of sensors 
and actuators: Ao = Ao (m 0 ,n). 

Actual air-fuel ratio A: Depends on several factors such as temperature de- 
pending enrichment during warm-up and correction for dynamic transients. 
At Diesel engines, lambda is always A > 1.3. 
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Figure 3.15 Piezoelectric injector map (IMPa = 10 bar) 



Battery voltage U b : It has an impact on the rise and fall times in electro- 
magnetically controlled injection nozzles. The effect can be compensated 
by adding a voltage dependant time correction A t(Ub). The compensated 
time is 

t inj + A t(U b ) . (3.36) 

Instead of direct measurement, the mass air flow m a can be estimated from intake 
manifold pressure p m or the throttle angle at at throttled SI engines. The air flow 
into the cylinders also depends on the dynamic pressure changes in the intake 
manifold. It is a function of 



= fo(Pm,Pm,n) , (3.37) 

where /o must be measured for all possible ?h 0 at stationary engine operation 
and corrected for dynamic pressure changes (Section 3.2.6). 

Figure 3.15 shows a map for the required injection time t,;„j of a high-pressure 
piezoelectric injector over the desired injection fuel mass mf and the injection 
pressure difference A p. 

3.2.5 Air Mass per Combustion Cycle 

The relative air supply X a at low engine speeds can be increased by acoustic 
resonances in the inlet pipes to each cylinder. These resonances are exited by the 
periodic opening and closing of the inlet valves. The geometry of the inlet pipes 
is designed for resonances at lower engine speeds. It is intended that a pressure 
maximum from the resonance occurs at the inlet valve at the time when it is 
opened. Hence, more air flows into the combustion chamber and increases the 
relative air supply A a and the effective work w e . Typical resonant frequencies are 
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m a 



4 
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Figure 3.16 Air mass m a obtained by integration of the air flow rh a over period tb — t a - 



which is simplified to: 



TO 



a 



T, 



CYL 



U \ (^1 ^ o ) u \ I 1 

m a (to) 7^2 m ay t i) 1 1 



— in a (tn) ^” + y 2 — +TO a (f„+i) ( 1 



(fl -ia) 2 \ 0 • /n \ 

y2 ) + 2 / y tra 0 (tj)- 

s ^ i=l 

(tn+1 - ib) 2 V 

7? A 



3.2.6 Intake Manifold Dynamics 

Figure 3.17 shows a cross-section of the intake manifold. The throttle angle 
controls the mass air flow ?h 0j j n into the manifold. Diesel engines are either 
unthrottled or very moderately throttled in some operating points in order to 
ensure a sufficient exhaust gas recirculation. Exhaust gas recirculation is not 
considered in the following model. 

The charge mass air flow out from the manifold into the cylinders m ac depends 
on the pressure level in the intake manifold p m (and the pressure in the cylinder 
p c ). To control the air- fuel ratio A correctly also in transients, the injected amount 
of fuel must be adapted to the mass air flow into the cylinder m ac rather than 
to the mass air flow ?h a y n into the intake manifold. 

The pressure oscillations shown in Figure 3.16 shall be neglected in the fol- 
lowing deduction (averaged model). A change in mass air flow m 0 results in a 
delayed change in manifold pressure p m . The according differential equation is 
derived from an energy equilibrium: The change of the internal energy of the air 
mass in the intake manifold is equal to the sum of in- and outgoing energy flows 
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plus the balance of energy changes of the gas due to the displacement work pV. 
By introducing the specific internal energy u = U/m and the specific enthalpy 
h = H/m the differential equation becomes: 

— (Tfl m U m ) — TTl a in'U'in 'Wac^air'U’ac T Pa7in PrrS^ac (3.42) 

at 

In this model, resonances within the manifold are neglected by assuming identical 
pressure p m at the throttle and the inlet pipes 1 . Additionally, the heat radiation 
from the engine is supposed to match the thermal heat required for the evapora- 
tion of the fuel. Therefore, no additional terms are added to Equation 3.42. The 
enthalpies at the inlet and outlet are equivalent to: 

Ala,inhin = Al a ^ n Ui n T PaYin (3.43) 

fll‘ac,airhac = ^ac^air^ac + PmVac (3.44) 

which can be inserted in Equation 3.42: 

= TTla^inhin 171 ac,air^ac (3.45) 

Using the specific heat coefficients c v = du/dd and c p = dh/dd as well as the air 

density p = m/V we get 

Pm b-m C-v^m H“ C-v^ mPm^m = 17la,inCpl^a ITlac^airCp'O rri (3.46) 



and after division by c v Vm 



q' - • q c p $ a 

Pm^m “ 1 “ Pm^m — T r 
Cv Vrr, 



171 a. in 






m„ 



(3.47) 



1 This relates to an averaged pressure model, which has proved to be sufficiently accurate 
for mass air flow transients. 
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Inserting the adiabatic exponent n = c p /c v and the ideal gas equation pV = 
mRiI, this yields the following equation for the pressure change: 



Pm — 



nRiIa 

V m 




n Wlac.air 
V a 



(3.48) 



It is difficult to measure the charge mass air flow from the manifold into the 
cylinder m a c,air ■ Because the dynamic response of m ac ,air is much faster than 
that of the manifold pressure p m , only the static behavior of rh ac ^ a i r shall be 
considered by a look-up table fi(n,p m ). The mass air flow m a c,air depends on 
the engine speed n and the manifold pressure p m . At stationary engine operation, 
where the pressure derivative is p m = 0, it can be determined by rh a ,in- 

Ift'ac,air — Plac,air — ,/l(?b Pm) — Pm — 0 (3.49) 

Pa 

The pressure change in the intake manifold is given by: 

Pm = -( rh a ,in ~ fi(n,Pm)) (3.50) 

T 

with the integration constant r: 



Vm 

nR/da 



(3.51) 



The look-up table can be measured on an engine test bed at stationary operating 
points, where the derivative of the averaged manifold pressure p m = 0. Under 
these conditions the incoming mass air flow m a) i n is equal to the values in the 
look-up table for rn* acair . Figure 3.18 shows the block diagram of the pressure 
model of the intake manifold. By integration of Equation 3.50 the average man- 
ifold pressure can be estimated in realtime. This can be utilized at engine idle 
speed control (Section 5.2). Another application is to estimate the mass air flow 
m a: i n from the measured manifold pressure p m and engine speed n: 



^o,m — 



T ■ p m + fl(n,p m ) 



(3.52) 
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The derivative p m can only be calculated with some lag. In some applications, 
the derivative of the throttle angle d* is therefore used as an additional variable. 
The time constant can be normalized at operating points p m , o and rh a o- 



T n 



Pm, 0 

— r 

m a , o 



(3.53) 



This normalized parameter r„ represents the integration time constant in seconds. 
The pressure change p m is then: 



d / Pm \ dfl a, in _ fl{n,Pm ) 

dt \Pm, o ) riiafi ?h Oj0 



(3.54) 



The time constant r„ depends mainly on the air flow fn a fi which can be seen 
in the following example: Assuming a manifold volume of V m = 4.25 1 at an 
ambient temperature of d a = 300° K and an adiabatic coefficient k = 1.4, the 
time constant r„ at minimum and maximum power can be calculated: 



• at maximum power, the manifold pressure is p m ,o = 1 b ar and the mass air 
flow shall be m a , o = 600 kg/h. This leads to a time constant of: 



T nj i = 21 ms 



(3.55) 



• at minimum power (e.g. idling), the pressure is p m , o = 0.35 bar and the 
mass air flow shall be m a , o = 6 kg/h which leads to: 

r n< 2 = 740 ms (3.56) 



It can be seen that the dynamic behavior of the intake manifold has an impact 
on engine dynamics especially at low power such as idling. 



3.2.7 Ignition Angle Control 

Correct ignition timing over the entire engine operating range is very impor- 
tant as it has a major impact on fuel consumption as well as on emission rates. 
Combustion within the cylinder can be divided into two phases: 

1. Inflammation delay (Time Proportional) 

In-cylinder pressure and temperature do not rise considerably within this 
time period. The inflammation delay time Tid depends on the temperature, 
pressure and air-fuel ratio. The delay can be convoluted into an equivalent 
crankshaft angle which increases with engine speed. 

2. Combustion (Angle Proportional) 

The equivalent crankshaft angle for the second phase is almost constant 
over the entire engine operating range. Induced by the piston movement, 
turbulences accelerate with engine speed, as well as the combustion process. 

If combustion starts too late because of retarded ignition angles, the emission of 
hydrocarbons HC will increase. High pressure amplitudes at advanced ignition 
angles increase the emission of NO x . NO x can be reduced by delaying the 
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ignition at the expense of a higher fuel consumption 2 . The following parameters 
are used to control the ignition angle: 

• Intake manifold pressure p m 

• Mass air flow m a 

• Engine speed n 

• Throttle angle a t 

• Air-fuel ratio A 

• Crankshaft angle acs and TDC signal of a reference cylinder 

• Ambient air temperature D a 

• Engine temperature D e 

• Battery voltage C4 

These values are the same as those needed for fuel control. The ignition angle on 
is dependent on many influences: 

• The ignition angle at is a function of engine load approximated by injection 
time ti n j ~ rh a /(n • A) (see Equation 3.33), and of engine speed n. This 
can be described by a map = f(U, n). The look-up table also covers the 
variation of inflammation depending on engine load and speed. Retarded 
ignition angles may be selected in order to compromise for reduced emis- 
sions and knocking. The ignition angle cq is determined for each engine 
operating point by test bed experiments. 

• Air- fuel ratio A, which determines the inflammation delay Tid- 

• Retarded ignition angle at high ambient air temperature D a to avoid knock- 
ing. A look-up table depending on and D a may be used. 

• Engine warm-up at low engine temperatures D e . A retarded ignition angle 
a-i retards the energy conversion process to a phase where the exhaust valves 
are already opened. The exhaust pipe and the catalytic converter are then 
heated very fast. 

• Engine speed stabilization at idling by advancing ignition angles at lower 
engine speeds, thus increasing torque. 

• Engine speed limitation by retarding ignition angles in conjunction with 
fuel cut off. 

• Retarded ignition angles during acceleration in order to avoid knocking. 

• Closed-loop knock control 

• The battery voltage C4 has an impact on the electrical ignition energy. 

Figure 3.19 shows an ignition angle map depending on engine speed and load. 
Angles on before Top Dead Center (TDC) are positive. 

2 Over all, the determination of the right ignition angle is a compromise between different 
objectives. 
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Figure 3.19 Ignition angle map 



3.2.8 Optimization of Engine Maps 

The fuel amount and ignition angle are the two most important parameters that 
influence the fuel consumption as well as the emission of pollutants. This has 
already been shown in Sections 3.1.2 and 3.1.6. There is a conflict between 
minimizing either fuel consumption or emissions. This is shown in Figure 3.20. 

If the ignition angle is chosen to minimize fuel consumption, the engine raw 
emission rates for NO x and HC will be fairly high. On the other hand, if the 
ignition angle is selected to minimize emissions, the fuel consumption will be 
higher. A compromise must consider fuel consumption and emission levels at 
all engine operating points. Emission levels can be very high at some particular 
operating points. There, the optimization must focus on the emissions. Other 
operating points show acceptable emission rates. At these points the optimization 
must focus on fuel consumption. 

Fuel consumption and emission levels are measured in special road driving 
cycles like the ECE-test or FTP-test. These tests specify the vehicle velocity 
over time. Translating vehicle to engine speeds, a test cycle is equivalent to a 
sequence of different engine operating points over time. Every operating point is 
defined by several control parameters including engine speed and load. 

The fuel consumption can be described by the volume V of combusted fuel 
over time. The minimization criteria is the integral over the test cycle. 

T 

J V ( t ) dt — > min 



V = 



o 



(3.57) 
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Figure 3.20 Fuel consumption and emission levels over ignition angle a;. 



The total fuel consumption V for a test cycle time T can also be obtained by a 
discrete summation over all visited engine operating points. 



N 

v = (3.58) 

i = 1 



An analysis of the test cycle shows that most operating points are visited several 
times. The individual time periods where the engine stays in the same operating 
point i can be summarized into a total time period t t . The fuel consumption 
over time Vi can then be minimized independently for each operating point. 
The resulting values of a,; and A,; are stored into look-up tables aii(ti n j,n) and 
A i{Unj,n) for every operating point. 

When optimizing fuel consumption, the maximum allowable emission levels 
are treated as optimization constraints. The maximum emission rates are fixed 
by laws which specify the maximum integral masses of the different pollutants 
generated during a test cycle. 
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N 

HC = ^HGia^X,) U <HC (3.59) 

i= 1 
N 

CO = ^2cO(cti,\i) ti <CO (3.60) 

i= 1 
N 

NO x = < iVO x (3.61) 

i=l 

The emission levels per time HC, CO, NO can be influenced by the values of 
tti and Xi at each operating point i. The emission limits are only given for the 
integral mass over the whole test cycle. It is therefore not obvious which cq 
and Xi values must be adopted at each operating point i. Such an optimization 
problem with constraints can be solved by using the Lagrange multiplication 
method [9]. The differences between actually achieved and acceptable emission 
levels are weighted by Lagrange factors L. Equation 3.58 and Equations 3.59 to 
3.61 are combined into a single criteria. 

W = V+L HC {HC-HC)+Lco{CO-CO)+L NOx (NO x -Nb x ) -> min (3.62) 

Now the cost function W must be minimized. For example, if all emission rates 
were at the acceptable limits (HC = HC, etc.), all terms except V for fuel 
consumption would disappear. This would minimize V as before. The value W 
can be divided into two parts: a constant part W 0 which is independent of the 
operating points i and a variable part influenced by cq and A *. 

N 

w = + 

i= 1 
N 

y^\L H cHC(ai, Xi)U + LcoCO(oti, X j)tj + L N o x NO x (ai, A*)t»] 

i- 1 

—LhcHC — LcoCO — Lno x N~O x 
N 

W = J2 Z ( a i’Xi)ti-W 0 

i—l 

where: 

Wo = LhcHC + L co CO + L NOx NO x = const (3.63) 

and: 

’Y*(®-i,X.i) V(cX’i,Xi') 

+ LhcHC (a,, A,) 

+LcoCO(a-i, Xi) 

~XL n o x NO x (oti, Aj) 
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Figure 3.21 Fuel consumption and NO x emissions over ignition angle cn and air- fuel 
ratio A at an operating point i with Lno x = 3. 



The value of W can be minimized by minimizing the function Z(oti,Xi) at each 
operating point: 

— * min (3.64) 

Figure 3.21 shows V (a*, A *) and NO x (ai, A*) for one operating point i depending 
on ignition angle on and air-fuel ratio \. The values were measured in an engine 
test bed run. The pair a,;, Ai for minimum Z(ai,Xi) is marked by a rectangle. 

It can be seen that fuel consumption V(ai,Xi) at minimum Z(ai,Xi) is 
slightly higher than the absolute minimum in order to compromise with emissions 
NO x (ai, Aj). The value W is a minimum for the entire test cycle if correct values 
for the Lagrange factors Lhc > Leo and Lnq x were selected. In practice, several 
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Figure 3.22 Look-up tables of optimum air-fuel ratio A and ignition angle a for an 
FTP test cycle. 



iterations with modified Lagrange factors must be passed to obtain a minimum 
W which also meets the legal emission constraints. For example, in a first iter- 
ation, low fuel consumption V is attained. However, emissions (HC,CO, NO x ) 
are still relatively close to or above the legal limits. The Lagrange weighting 
factors L are then increased for the next iteration. Practical experience shows 
that only a few iterations are needed to minimize the function Z at a small 
number of representative operating points i. Values for other operating points 
can be obtained by interpolation. Figure 3.22 shows resulting look-up tables for 
the ignition angle a and the air-fuel ratio A. In engines at stoichiometric A = 1 
operation, only the ignition angle on is calculated. The optimization approach 
can also be extended to other variables such as exhaust gas recirculation rate. 




4 Diesel Engine Modeling 



In recent years engine modeling has become more and more important in the 
development process. In an early development stage new strategies can be exam- 
ined without expensive test bench measurements. Further on, control strategies 
can be designed and optimized. Thereby new requirements for the models arise. 

• Multiple direct fuel injection in diesel engines 

Direct fuel injection was a breakthrough for diesel engines. By means of 
multiple injections in common rail systems, the evaporation process and 
subsequently the combustion process can be directly influenced. 

• Stratified charge in SI engines 

As explained in Section 5.1, SI engines work with a homogeneous air-fuel 
mixture at stoichiometric air-fuel ratios because the catalytic converter 
shows extremely good conversion rates under such conditions. At part 
load, directly injected SI engines work with lean air-fuel mixtures. In or- 
der to get a safe combustion a stratified charge is assembled around the 
spark plug. In such engines the charge exchange, turbulences inside the 
combustion chamber and the injection process must be modeled. 

• Exhaust gas calculation 

One of the major problems is to model the exhaust gas generation fairly 
accurate. The maximum allowable emissions are determined by law. In 
order to meet such regulations two strategies are pursued in parallel. One 
is to minimize the raw emissions of noxious gases from the combustion 
process. The other is to clean the exhaust gases by means of catalytic after 
treatment systems. Both strategies need sophisticated control systems and 
dynamic models of the chemical processes involved. 





76 



4. DIESEL ENGINE MODELING 



Intake Stroke Compression Stroke Power Stroke 




Exhaust Stroke 
exhaust 

pi-u-ds 



/ 



0 



Figure 4.1 Four stroke cycle diesel engine. 



Because of the rising importance of diesel engines, models for flexible injection, 
energy conversion and evaporation processes will be presented in this chapter. In 
addition a model for soot accruement is introduced. 



4.1 Four Stroke Cycle Diesel Engine 

The movement of the piston up and down inside the cylinder is called its stroke. 
The crankshaft is connected to the piston by the connecting rod. By this means 
the linear motion of the piston is converted into rotational movement. The crank 
shaft connects all the pistons and the engine torque. Most diesel engines use a 
four stroke cycle because the piston has to travel up and down twice to complete 
a full power cycle turning the crank shaft through 720 degrees. As listed below 
the four strokes of the diesel cycle are illustrated in Figure 4.1. 

• Intake stroke 

The intake stroke begins at top dead center. As the piston moves down, 
the intake valve opens. The downward movement of the piston draws fresh 
air through the intake valve into the combustion chamber. 

• Compression stroke 

The compression stroke begins with the piston at bottom dead center. The 
piston is rising up to compress the in-cylinder charge. Since both the intake 
and exhaust valves are closed, the in-cylinder gas is compressed to a fraction 
of its original volume heating it up. 

• Power stroke 

The next stroke is called power stroke. At the start of this second downward 
movement of the piston or at the end of the compression stroke fuel is 
sprayed into the cylinder. On contact with the hot in-cylinder charge the 
fuel inflames and the resulting combustion drives the piston down. 
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• Exhaust stroke 

The last stroke is called exhaust stroke. Here the exhaust gases are pushed 
out of the combustion chamber through the outlet valve. 

In the following equations for the four stroke diesel engine model, the crank- 
shaft angle acs is taken as an independent variable instead of time t. The 
relationship between these variables is 

acs = u> ■ t = 2n n-t , (4-1) 

dacs = u> • dt , (4-2) 

with n as the engine speed. 



4.2 Charge Exchange 

During charge exchange the exhaust gases inside the combustion chamber are 
replaced with fresh air. The volumetric mass flow through a restrictor with 
the cross sectional area A e ff into and out of the combustion chamber can be 
calculated according to [82] as 
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tHiov, A eff ^ ^ 


K ~ 1 


\vo) 


\PoJ 



— Pflow • 7l e // • — q- • 4/ . (4-3) 

vR ' 

The outflow 4/ depends on the pressure po in front of and on the pressure p\ 
behind the throttle as well as on the adiabatic exponent k. In practice, the 
effective cross sectional area is smaller than the geometric cross sectional area. 
This fact is considered in the model with the flow ratio coefficient Pfi ow • The 
temperature before the throttle is d 0 . 

The geometric cross sectional area A e f f depends on the valve movement and 
may be described as 

A e ff = n a ■ 7r • h v • cos(cr) • ( d v + h v ■ sin(cr) • cos(cr)) , (4.4) 

depending on the number of input or output valves n 0 , the valve stroke h u , the 
valve seat angle a and the inner valve seat diameter d v . Figure 4.2(a) shows the 
geometric parameters at the valve seat. The valve stroke h v during inlet and 
outlet can be approximated with sine curves, see Figure 4.2(b). 

4.2.1 Flow into Exhaust Pipes 

After the combustion, the exhaust gases are pumped into the exhaust pipes. 
At the beginning of the outlet phase the in-cylinder exhaust gas approximately 
consists of a homogeneous mixture of air and burnt fuel. The air mass m a and 
burnt fuel mass mfj )urn t have to be adopted after completion of the combustion 
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(a) Geometric parameters (b) Valve stroke 

Figure 4.2 Cross sectional area at the valves. 



process. The flow of exhaust gases through the outlet valves is depending 

on crankshaft angle instead of time. It is calculated according to Equation 4.3. 




(4.5) 



In this, Pout is the pressure in the outlet pipe, p the in-cylinder pressure and D 
the in-cylinder temperature. Integration of Equation 4.5 over the outlet phase 
yields m ou t . Some remaining gases m rem will remain in the combustion chamber 
after closing of the outlet valves. 

m rern : m a T .burnt ^71 out (d-6) 

This remaining gas consists of air m re . m ^ a i r (see Equation 4.14) and burnt fuel 
mrem, burnt (see Equation 4.15). Since the exhaust gas is assumed to be a homo- 
geneous mixture, the air-fuel ratio of the remaining gases is equal to that of the 
gas mixture at the end of the previous combustion phase. 

4.2.2 Flow into Combustion Chamber 

The mass flow drriac into the combustion chamber from the intake manifold is 
aacs 

calculated with Equation 4.3 in dependency of the crankshaft angle acs and the 
engine speed n as 




(4.7) 



In this, p m is the intake manifold pressure and p the in-cylinder pressure. Fur- 
thermore, the gas composition must be considered in the charge process, such as 
in-cylinder remaining gases m rem and recirculated exhaust gases itiegr (exhaust 
gas recirculation, EGR). By means of a given EGR-rate Xegr, the charge air 
mass from the manifold is 



m ac , a i r — m ac • (1 Xegr ) 



(4.8) 
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and the mass of the recirculated gases is 

mEGR = m-ac ■ X EGR ■ (4.9) 



The charge mass of fresh air per stroke m ac is calculated by integrating Equation 
4.7 over the inlet phase. Since the recirculated gas is composed of air and burnt 
fuel at the same proportion as the remaining gases at the end of the preceding 
outlet phase, the amount of recirculated air is 



TnEGR,air = 



TO. 



rem.air 



■ X 



m r 



EGR ■ m a 



(4.10) 



The resulting air mass m a inside the combustion chamber for the next combustion 
phase is the sum of the air mass from the manifold rn ac ,ain the recirculated air 
mass to EGR,air and the air portion of the in-cylinder remaining gases m r em,air 
from the outlet phase 



TOa = m ac ■ (1 — XeGR ) + RT'EGR.air + TYlrem,air • (4.11) 

The amount of recirculated burnt fuel m EGR , burnt, is the complement to equation 
4.10 

fUEGR, burnt = X EGR ■ TO oc - m EGR ,air = X EGR • TO ac • — rem ’ burnt _ (4.12) 

"Irem 

The resulting charge mass of burnt fuel burnt a t the beginning of the combus- 
tion process is the sum of recirculated burnt fuel m EGR , burnt and the portion of 
the in-cylinder remaining gases m rem ^urnt from the previous outlet phase. 

TOy , burnt = RT'EGR, burnt A RT'rem, burnt • (4-13) 

The injected fuel is not yet burnt at the time of injection. Figure 4.3 shows all 
masses determined during charge exchange. 



4.3 Air-fuel Ratio 

The stoichiometry ratio, L st refers to an air-fuel ratio in which all combustible 
materials are used with no deficiencies or excesses. The air-fuel ratio A c is the 
current ratio of air and burnt fuel masses divided by the stoichiometric constant. 
Local differences of the air-fuel ratio are neglected in this model. The gas is 
assumed to be homogeneous. 



4.3.1 Exhaust Stroke 



In this section of the diesel cycle, the air-fuel ratio A c of the in-cylinder exhaust 
gases shall be identical to the air-fuel ratio at the end of the previous combustion 
process (section 4.3.3). Since we assume the exhaust gases to be homogeneous, A c 
remains constant during the exhaust stroke. The resulting in-cylinder remaining 
air mass ?n rem , a ir and burnt fuel mass m rern ,, burnt, are now derived as 



1 

"Irem, burnt — T , ' T ' ^rem 

A c ‘ Lst + -L 
* Lst 

Til rent, air T j . T ' Tll rern . 

A c ‘ L s t + -L 



(4.14) 

(4.15) 
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Figure 4.3 Overview of all calculated masses during charge exchange. 



4.3.2 Intake Stroke 



Before fuel injection the combustion chamber is filled with fresh air and recircu- 
lated gases. Since the combustion chamber is charged with fresh air, the air-fuel 
ratio A c is increasing. The resulting air-fuel ratio is 



Ac = 



burnt * I'st 



■ (1 — XeGr) + rn.EGR,air + m r 

il^EGR.burnt T ^Vem, burnt) ' -C 



(4.16) 



m re m, a ir and m re m, burnt are the remaining gas masses of the previous exhaust 
stroke (section 4.3.1). The change of the air-fuel ratio A c during the intake stroke 
is 



d\ c 1 ( dm a dm.f bur nt\ 

7 — o r 1 j * TTl f , burnt m a * , I 

dacs rn /,i burnt L st \da C s dacs ) 



(4.17) 



4.3.3 Compression and Combustion 



The air-fuel ratio A c is changing during combustion as well. Changes of the gas 
masses inside the combustion chamber result from the burning of injected fuel 
which causes the air-fuel ratio to decrease. The current ratio A c between air and 
burnt fuel is calculated as 



y _ ITT-a _ ITT-ac ' (1 ~ X E Gr) + 1REGR,air + TTlrem,air ^ 

ITlf, burnt ’ Lst {iTl EG R, burnt H"~ ^rem, burnt ^in j , burnt') ' Lst 
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m a represents the mass of fresh air and the mass of air contained in the com- 
bustion products inside the cylinder. After closing the inlet valves m a remains 
constant during the compression and the power stroke. The amount of burnt fuel 
to /, burnt > consisting of the amount of burnt fuel in the remaining gases m rern ^ U rnt , 
the recirculated burnt fuel mass rriEGR, burnt and the amount of injected and al- 
ready burnt fuel 'rrii n j. burnt changes during the current combustion process. The 
change of the air-fuel ratio A c in dependence of the crankshaft angle is 



d\ c 

dotes 



A c CLTTlin j , burnt 

^ / , burnt dacs 

Ae 

RT'EGR, burnt T Rblrem, burnt T R^inj , burnt 



dm 



in j , burnt 



dacs 



. (4.19) 



The air-fuel ratio at the end of combustion is taken in Equation (4.14) and (4.15). 



4.4 Mass Balance 



A gas mass increment dm. inside the cylinder is the sum of the charge mass 
increment dm ac and the increment of the injected fuel that is already burnt 
dmi n j, burnt minus the outlet mass increment dm ou t- Blow-by effects at the valves 
and at the piston may be neglected for modern diesel engines. 

dm dm ac ^ d J m- tn jbumt dm^t ^ j 

dacs dacs dacs dacs 

During the intake stroke, the gas mass derivate is 



dm. dm ac 
dacs dacs 



(4.21) 



During the compression and power stroke, the valves are closed and the mass 
increment dm is the combusted fuel mass increment dmi n j^ U mt 



dm dmiuj^burnt 

dacs dacs 



(4.22) 



During the exhaust stroke the mass derivative depends on the flow of exhaust 
gases out of the cylinder. 

dm = dmout ^ 

dacs dacs 

The liquid fuel inside the combustion chamber is ignored in this balance since 
only gases are relevant in the combustion process. 



4.5 Fuel Injection 

The energy conversion is determined by the injection process, i.e. how much 
fuel is injected over time. The course of injection is influenced by the injection 
pressure, the number of injections per combustion cycle, the duration of each 
injection, the geometry of the injector nozzle and the in-cylinder pressure. 
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\ 



fuel supply 



Figure 4.4 Schematic structure of an injector. 



Figure 4.4 shows the schematic structure of an injector. The fuel inside the 
injector is under the same high pressure as inside the common rail. Because of the 
needle movement, the nozzles open and fuel flows into the cylinder according to 
the pressure gradient between common rail and combustion chamber. The flow 
of the injected fuel mass da J s in dependency of the crankshaft angle acs can be 
calculated with the Bernoulli Equation for incompressible substances [98, 121] as 



drrif 

dacs 



/J ' mJ ' A n ' J2 ■ pf ■ (j) rai i - p ) 
•7 T ■ n V 



(4.24) 



The injected fuel mass over rotation of the crankshaft depends on the effec- 
tive cross sectional area of the nozzles A n , the fuel density pf and the pressure 
difference between the injection pressure p ra ii and the in-cylinder pressure p. 
Furthermore, the effective cross sectional area A n of the nozzles depends on the 
needle movement. The factor p,i n j is required for fitting calculated fuel masses 
to actually measured fuel masses in each operating point. Because the needle 
movement is often not known, a model without needle movement is employed. 
According to [121] the fuel flow into the combustion chamber for direct injected 
diesel engines can be assumed as parabolic over rotation of the crankshaft. 

The parabola 



dm f (a C s) 
dacs 



a ■ a 2 cs + b ■ acs + c 



(4.25) 



has to be open downward, i.e. a < 0. The maximum point (acs, max, is 

calculated as 



acs, max — 0 and 

2 • a 

The injection flow can be written 



dmf tmax _ 4 ac - b 2 
dacs 4a 



dm f (a C s) 
dacs 



a ■ a 2 cs + b ■ a C s 



0 < acs < Cncs,inj j 



(4.26) 



(4.27) 



if the injection starts at acs = 0- The end of injection occurs at angle acs,inj- 
The maximum fuel flow results from Equation 4.24 at maximum cross sectional 
area A n ^ max 



dm max 

dacs 



fJ'inj 

2 • n ■ n 



■ A 



n,max 



\J%- Pf ■ (Prail ~P)- 



(4.28) 
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The crankshaft angle of maximum fuel flow is 

OtCS, ini 

&C S,max — q 



(4.29) 



The injection process is finished at acs = oics,inj ■ The parameter a is calculated 
at the end of the injection, where is again zero in Equation 4.27 

a = ~ b . (4.30) 

&CS,inj 



Using the maximum flow in Equation 4.28, angle acs, max and Equations 

4.27 and 4.30, the parameters of the parabolic equation are found as 

a = ' O /imJ ' A n,max ' J % ' Pf ' ( Prail ~ P) (4-31) 

a CS,inj 2 • 7T • n v 

b = • /i ” iJ • A n max ■ J 2 ■ pf ■ ( Prail -p). (4.32) 

acs,inj 2-7 T-n V 

Integrating Equation 4.24 over the injection angle segment acs,inj yields the 
total injected fuel mass irif for one combustion cycle 



_ 4 Pini , 

m f = Cr^' A n 

&CS,inj 



&CS,inj 



2 • Pf ' ( Prail ~ P ) ' 



a cs 

OtCS,inj 



acs - ~ — ) dacs- 
(4.33) 

The integration is done stepwise, generating a fuel package m/ : ; for each angle 
segment A acs within limits 0 < A acs ■ i < acs, ini- 



4 * A n rnax Finn , 

m ‘ J = OCSM ■2^'P'l>r(Pr««-p)- 



A otes-i 



A QCS'(>-1) 



acs - 



a cs 



acs,i 



dacs- 



(4.34) 

If the entire injected fuel mass m/ is given, the injection angle segment acs, ini 
is calculated as 



3-7 T-n 



acS, ini ~ m f ‘ 



* -An, max ■ V 2 ’ Pf ’ ( Prail ~ P ) 



(4.35) 



To parameterize the fuel injection model, a map for pi n j is determined for differ- 
ent operation points. Figure 4.5 shows a comparison of measured and modeled 
data of the injection time t-i n j = ^^'n 1 for several pre- and main-injections. 



4.6 Fuel Evaporation 

The piezo electric injection technique allows to shape the course of injection and 
the progress of fuel evaporation. Many different injection curves are possible, 
see Figure 4.6. The pre-injection helps to reduce the typical hard diesel engine 
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Figure 4.5 Comparison of measured and calculated injection durations for pre- and 
main-injection. 



injection curve 




combustion noise. The pressure gradient can be flattened with a small fuel in- 
jection at the beginning of the combustion cycle. The main injection and the 
shape of the main injection are used to optimize the power output of the engine. 
Finally, post-injections are useful for exhaust gas after treatment systems, e. g. 
to heat up the exhaust system and to burn soot in a catalytic converter inside 
the exhaust pipe. 

For this reason the simulation model must be able to emulate this behaviour. 
Simple energy conversion models like the Vibe- or the Double-Vibe-functions 
[86, 98, 125] are not sufficient. Rather a phenomenological approach according 
to Constien [20] is employed in this book to model the energy conversion course. 

During injection, the fuel sprays into the combustion chamber and evaporates 
depending on the injection pressure, the fuel density, injector nozzle geometry 
and the in-cylinder pressure. The air-fuel ratio is inhomogeneous inside the com- 
bustion chamber. Zones with lean and rich mixtures develop during injection. 

The injected fuel evaporates into drops and droplets. Around the drops, a 
gaseous environment develops. The temperature is increasing from the center of 
each drop to its surface where it reaches the combustion chamber temperature. 
In addition to pressure, the self inflammation of evaporated fuel depends on the 
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(a) Theoretical injection curve 





■ liquid 

□ gaseous 

□ burnt 

(d) Legend 



Figure 4.7 Schematic sequence of the evaporation process. 



temperature and the air-fuel ratio around the droplet. When the combustion 
starts, the droplets heat up and the oxygen concentration is decreasing. Around 
the droplets lambda values between zero and infinity occur. 

According to [20, 122] the number of droplets, their surface and the amount of 
evaporated fuel inside the combustion chamber can be estimated. This informa- 
tion is needed to calculate the energy conversion process. Figure 4.7 shows how 
the evaporation process is modeled. An exemplary injection curve with pre- and 
main-injection is depicted in Figure 4.7(a). The fuel is assumed to be injected 
into the combustion chamber in discrete subsequent fuel mass packages m.f^. 
For every fuel package the portion of liquid, gaseous and burnt fuel develops over 
time. With the crankshaft angle as the independent variable, the condition of 
each fuel mass package m/y is therefore calculated over angle segments A acs • j- 

Figure 4.7(b) shows a snapshot of the fuel masses in the combustion chamber 
at the momentary crankshaft angle acs, i- The perpendicular lines separate the 
different fuel mass packages rrifj . Within each package, the three fuel portions 
are shown. Different gray tones mark the liquid, gaseous and burnt fractions of 
the injected fuel. 

At the momentary crank angle acs , 2 almost all fuel has been injected, Figure 
4.7(c). It can be seen that most of the fuel is already burnt, especially in the 
packages of the pre-injection. The fuel package singled out at acs . l is almost 
completely burnt at acs, 2 - Only a little gaseous fuel remains. 

The number of drops and the surface area of the drops shall now be calcu- 
lated. Because it is impossible to analytically describe the multitude of all drop 
diameters, an average drop diameter d 32 in dependence of injection pressure and 
fuel package mass is adopted. According to the Sauter approach [20], it is cal- 
culated as the ratio of the entire drop volume divided by the entire drop surface 
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area 



^32 — 



E d 3 dr -H r 
Ed 2 dr -H r 



(4.36) 



The number of drops with the same diameter d dr is H r . Since all d dr and H r are 
unknown, the Sauter diameter ^32 may alternatively be calculated according to 
the Varde/Popa/Varde [45] approach 



d 32 = 16.58 (Re ■ We ) -0 ' 28 



(4.37) 



where Re is the Reynolds number and We the Weber number of the flow. The 
Reynolds number Re is calculated depending on the fuel flow velocity at the 
injector nozzle v n , the nozzle diameter d„ and the kinematic fuel viscosity Pf 

Re = . (4.38) 



The Weber number We is calculated with the density of the combustion chamber 
charge p charge and the fuel surface tension 07 as 



We = 



Pcharge 






(4.39) 



To get the fuel flow velocity v n at the nozzle, the pressure difference A p be- 
tween the injector and the combustion chamber as well as the fuel density pf are 
required 



Vn = 




(4.40) 



Inserting Equations 4.38, 4.39 and 4.40 into Equation 4.37, the Sauter diameter 
d 3 2 can be written as 



d 32 = 12.392 



* 44 wr • (-/ • ^/)°- 28 

Ap °- 42 • P° ch 2 a 8 rge 



(4.41) 



The number of drops N dd in a fuel package i injected at the discrete crankshaft 
angle A acs • * is derived from the injected liquid fuel mass to/, * 



Naa = 



mfA 

f ' ^32 ' Pf 



(4.42) 



TO /4 is derived from the angle-discrete integration of Equation 4.34. The number 
of drops N dd is calculated once for each fuel package. The portions of liquid, 
gaseous and burnt fuel in the package i change over the discrete crankshaft angle 
progression A acs ■ j ■ Also the drop diameter d 3 2 ,ij changes over subsequent 
angles Ao^s • j- The drop diameter in package i is recalculated iteratively as 



^32 ,ij 




m f,jj 
■ Nd,i ■ Pf 



(4.43) 
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At the crank angle A acs ■ j = A acs • i, the liquid fuel mass has the initial value 
m,f ij = rn f -i (Equation 4.34). The surface area of all drops in package i at the 
angle A acs • j is 

= Nd,i ■ 7T • ^32 ,ij ■ (4.44) 

The evaporated fuel mass of package i at the angle A acs ■ j is 

A m fiev ,ij = C di ff ■ A fi ij ■ p mp ■ Aa / cs , (4.45) 

n ■ Cl32,ij 

where p is the combustion chamber pressure, n the engine speed and Cdiff the 
diffusion constant which is determined in Section 4.8. By means of exponent m p 
the dependance of the evaporation process on the in-cylinder pressure can be 
fitted. Variable A acs is the discrete crankshaft angle step size. 

The liquid fuel portion of fuel package i at the next step j + 1 is calculated 
by subtracting the evaporated fuel Am/. elV y of Equation 4.45 from the liquid 
portion m f tl j of the last step j 

mf t i j+ 1 = m,f t ij - A m/ >ev ,»j • (4.46) 

The drop diameter <1^.2, i j+ 1 an d th e surface area Af i J+1 are calculated from 
Ecpiation 4.43 and 4.44 when increasing j to j + 1. 

In order to describe the transformation of evaporated into burnt fuel the 
inflammation delay time t^j is determined at each step j. The inflammation 
delay time Tid.j nray be approximated by the empirical approach [20] depending 
on the average in-cylinder pressure p.j and the average temperature i)-j between 
A acs • i an d the momentary angle segment A acs ■ j 

2100 K 

ndj = 2, 1 • pL ’ ■ e . (4.47) 

After the inflammation delay Tidj, the evaporated fuel mass m.f :eVt ij of fuel 
package i at the beginning of iteration step j and the evaporated fuel mass 
A m.f tev ,ij during step j are assumed to burn completely into 

A TTl f 'burnt ,ij — m f,ev,ij -)- Nnzf eV’ij • (4.48) 

This occurs at the angle segment A acs'j > Acres ' * + 27 t « ■ Tid.j ■ Before the next 
iteration step j + 1, mf teVt i j+i is set to zero, because new fuel will evaporate at 
the step j + 1. If the inflammation delay is not yet reached the overall evaporated 
fuel mass j+i of package i is updated at the beginning of the next iteration 

step j + 1 

TYlf,ev,i j+1 = 'lTlf,ev,ij “I” A7Tlf^ e v,ij • (4.49) 

The entire burnt fuel mass mf^umt,j at iteration step j is the sum of all burnt 
portions of the different injection packets i 

ITl burnt, j — E ^Tfl f burnt ^ij (4.50) 

i 

and the mass of the injected and already burnt fuel rrii n j y burnt is the sum of the 
already burnt fuel masses 

'Wlinj, burnt — E ^ f, burnt, j • (4.51) 

j 
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Figure 4.8 Charge and combustion model. 



4.7 Cylinder Dynamics 

The aim of modeling the combustion process is to get information about the in- 
cylinder conditions for different operating points of an engine. The development 
of engine control strategies and of the calibration procedure may be simplified 
with such models. Depending on the requirements, different models can be used. 
One application of the zero-dimensional model are for hardware-in-the-loop con- 
figurations. A more complex two-dimensional model can be found in [122]. 

Both models base on the ideal thermodynamic gas Equation 2.1. For the 
calculation of the combustion process itself, the thermodynamic balance 

dwt + dQ + dm r ■ ( h r + e r ) 

r 

= dU + dE ext (4.52) 

is regarded. Wt is the technical work, Q the external heat, dm r are different 
infinitesimal masses crossing the system border (combustion chamber) with the 
enthalpy h r and specific external energy e r . 

The right side of Equation 4.52 represents the energy stored inside the system, 
where U is the internal energy and E ext the external energy e.g. kinetic or 
potential energy of the system. 

Both approaches employ the phenomenological evaporation model in Section 
4.6 and the fuel- inject ion model in Section 4.5 for flexible energy conversion rates. 

4.7.1 Zero-Dimensional Modeling 

In zero-dinrensional modeling, the combustion chamber is regarded as a closed 
system with a single reaction zone. The gases inside are assumed to be ideally 
mixed. Local differences of the calculated variables, flows or eddies are neglected 
in the zero-dinrensional model. Figure 4.8 gives an overview of the combustion 
chamber with some of the state variables. 

In the following, the thermodynamic equations are formulated in differential 
notation. In order to calculate the energy balance in all four strokes of an engine 
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cycle, the energy conversion during combustion including heat transfer, charge 
exchange and volume change is integrated. 



4.7.2 Thermodynamic Equations 

The ideal gas equation 2.1 is differentiated with respect to the crankshaft angle 

«cs 



p- 



dV 

dacs 



- V • 



dp 

dacs 



= m ■ R ■ 



dd 

da cs 



+ m • d • 



dR 

da C s 



+ R-d • 



dm 

da C s 



(4.53) 



According to [98] and [111], a variation of the gas constant R influences the 

thermodynamic state equations only for temperatures above 1800A' or for air- 

fuel ratios below 1.2. In diesel engines the combustion temperatures are usually 

below 1800iv and the air-fuel ratio is above 1.4. Hence, the gas constant R is 

considered to be constant and , dR = 0. 

dacs 



4.7.3 Energy Balance 

According to Equation 4.52 the internal energy U inside the cylinder changes if 
there are changes of the external heat dQ , the technical work dw t or the enthalpy 
dH crossing the system border. Because of its minimal impact, the influence of 
the external energy dE ext is neglected. 



dw t 

dacs 



dQ w 

dacs 



- h 



out 



dm, out 
dacs 



+ h ac ■ 



dm ac 

dacs 



dQ comb 



dacs 



dU 

dacs 



(4.54) 



is the volumetric work of the piston. The heat loss through the wall is 

laCs ' ^ out ' d da°cs ent halpy change of the mass flow out of the combustion 

chamber. Equivalently the enthalpy change of the mass flow through the inlet 
channel into the cylinder is h ac ■ dacs ■ is the released energy during 

combustion. All this sums up to the change of the internal energy of the gas 

dU 

dacs ' 



4.7.4 Volumetric Work 

During the combustion, the combustion chamber is a closed system. The volu- 
metric work drives the piston. Equation 2.4 relates volumetric changes to work. 
To calculate the volumetric work, the piston position is necessary. In Equation 

3.4 the piston position s in dependency of the crankshaft angle acs is given. 
At Top Dead Center (TDC), the volume is minimal, the so called rest volume 
Vtdc- The current volume is calculated from the piston surface A piston , the total 
displacement volume Vd and the piston position s(acs ) 



V(acs) 



Vtdc + s(acs) • A pist on 



Vd 

£ — 1 



+ s(acs) ■ 



Vd 

2 ■ r ' 



(4.55) 




